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Abstract 

In the previous study[H [2], we formulate a matrix model renormalization group based on 
the fuzzy spherical harmonics with which a notion of high/low energy can be attributed 
to matrix elements, and show that it exhibits locality and various similarity to the usual 
Wilsonian renormalization group of quantum held theory. In this work, we continue the 
renormalization group analysis of a matrix model with emphasis on nonlocal interactions 
where the fields on antipodal points are coupled. They are indeed generated in the renor¬ 
malization group procedure and are tightly related to the noncommutative nature of the 
geometry. We aim at formulating renormalization group equations including such nonlocal 
interactions and finding existence of nontrivial held theory with antipodal interactions on 
the fuzzy sphere. We find several nontrivial hxed points and calculate the scaling dimen¬ 
sions associated with them. We also consider the noncommutative plane limit and then 
no consistent hxed point is found. This contrast between the fuzzy sphere limit and the 
noncommutative plane limit would be manifestation in our formalism of the claim given 
by Chu, Madore and Steinacker that the former does not have UV/IR mixing, while the 
latter does. 
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1 Introduction 

Matrix models have been considered as constructive definition of quantum gravity. Rela¬ 
tively simple one- and two-hermitian matrix models have been solved exactly and provide 
nonperturbative formulation of two-dimensional Euclidean quantum gravity coupled to 
c < 1 conformal matters M- Yang-Mills types of matrix models with fermionic symme¬ 
try have been proposed as nonperturbative formulation of string/M theory [5]. Each of 
them is given as a lower dimensional reduction of maximally supersymmetric Yang-Mills 
theory and have not yet been solved exactly as in the case of two-dimensional quantum 
gravity. 

The large-A renormalization group (RG) method proposed by Brezin and Zinn-Justin 
[6] is an analytic approach to investigate critical behavior of such a matrix model in the 
large-A limit. For simple one- or two-hermitian matrix models where exact results are 
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available, it has been shown that large-iV RG analysis captures critical behavior of the 
models, first qualitatively j.6] and later even quantitatively with the help of Schwinger- 
Dyson equations [7J. 

Though the original motivation of this large -N RG method is to analyze models that 
we cannot solve exactly, and to find a way to overcome the so-called c = 1 barrier, the 
actual application has not been explored much so farQ In the past decades, the focus 
of the study on large -N matrix models as quantum gravity has been shifted to reduced 
supersymmetric Yang-Mills type matrix models as mentioned above. There have been lots 
of study on these models (see e.g. [ID]), including very impressive numerical studies (e.g. 
HU), but it is of course very plausible if we can develop an analytic approach to dissect 
such models. One of the characteristics of these “new” matrix models is that the matrix 
elements have direct physical interpretation; for example, they are often interpreted as 
positions of D-branes and open strings connecting them. On the other hand, in the case of 
matrix models for Euclidean D < 2 quantum gravity, the matrices are technical tools to 
sum up random surfaces with spin degrees of freedom that give rise to conformal matters. 
Thus, the large-iV RG method proposed by Brezin and Zinn-Justin somehow inherits 
this nature; namely, there are no criteria on which matrix elements are to be integrated 
out. On the other hand, in the modern interpretation, we may want to attribute some 
physical meaning to matrix elements, and formulate large -N RG based on this. With 
this motivation, together with D. Tomino, we have formulated a large -N RG by using 
the fuzzy spherical harmonics [12J. The fuzzy spherical harmonics are an analogue of 
the standard spherical harmonics, and they span a basis on which general hermitian 
matrices can be expanded. Thus, there appears a clear notion of high/low energy modes 
for expansion coefficients. We can in this way formulate the large -N RG, and have shown 
that they enjoy several nice properties such as locality and derivative expansion of double 
trace terms. We also find Gaussian and non-Gaussian fixed points, and discuss their 
properties. However, we have also faced a difficulty that originates in a noncommutative 
nature of the geometry; it is embodied as emergence of new nonlocal interactions between 
fields on antipodal points on the fuzzy sphere. In the previous work, we carry out the RG 
analysis by dropping these antipodal interactions. In this paper, we come back to this 
issue and present more complete treatment of RG analysis. 

In the following subsection, we start our discussion with what we have done in the 
previous paper and what we want to do in this paper. 

1.1 Large -N renormalization group based on fuzzy spherical 
harmonics 

We first review the basic formulation of the large -N RG based on the fuzzy spherical har¬ 
monics!]] We then describe its intriguing aspects that have been observed in the previous 
study [1]. 

*One may see mm for some trials. 

tfiere only basics which are necessary to explain the problems are presented. More detailed explanation 
is provided later. 
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We start with a matrix model that describes 0 4 scalar field theory on a fuzzy S 2 of 
radius J3 


S N = P ftr K ^ + ^A , (1.1) 

where N denotes the size of the hermitian matrix 0, L* is the generator of the SU(2) 
in the spin L representation with N = 2L + 1, and m 2 N and /% are the mass (squared) 
and the coupling constant respectively^ 0 is expanded by use of N x N fuzzy spherical 
harmonics matrices, 


2 l i 

0 = AmTlm- ( 1 - 2 ) 

1=0 m=—l 

The maximum of the angular momentum 2L is related to the size of the matrix N as 
N = 2L + 1. The expansion coefficients 0; m are considered to be dynamical degrees of 
freedom in field theory of 0, and thus there exists a clear notion of “high energy” modes. 
The large-iV RG transformation is defined by integrating over 2L +1 number of maximum 
angular momentum modes, 02 Lm with — L < m < L. This can be carried out by usual 
perturbation theory, and it generates tr0 2 and tr0 4 terms (of the reduced size). Upon a 
suitable change of the matrix basis, these changes can be absorbed into the mass and the 
coupling constant. This is an analogue of Wilsonian RG flow, and we can find Gaussian 
as well as non-Gaussian fixed points in the previous study [1]. In fact, this procedure can 
also be regarded as the Wilsonian RG of a scalar field theory on the fuzzy sphere. We 
also examined linearized RG transformation around the fixed points and calculated the 
scaling dimensions of operators at the fixed points as its eigenvalues. Various large-iV 
limits with different noncommutativity respected are also argued. 

The perturbative calculation of the RG transformation also generates various terms 
that are not present in the starting action (11.11) : up to C>((fi;®) 2 ) we considered, it includes 
trjv0 6 term and also double trace terms. It turns out that in the low energy regime 
l L these terms are either negligible or can be rewritten as a series of single trace 
operators with derivatives (commutators with L0s) which are suppressed in 1/N. Namely, 
it has been observed that the corrections are controlled as in the case of usual Wilsonian 
RG analysis. However, it is found that there appears another type of correction terms 
that does not present in the usual Wilsonian RG. Let us consider the following simplest 
example, 


2L-1 i 

tr N (00 out 00 out ) oc ^2 E (-A l <t>imtr N (ct>T lm ) 

1=0 m=—l 

Un Appendix [3 we summarize the basic properties of the fuzzy sphere and also construction of the 
matrix model for scalar field theory on it. 

Tn PQ, k'A was called gjy. In this paper, the symbol g n is reserved for other combinations of the 
coupling constants. In later sections, stands for the coupling constant for the same operator tr at0 4 . 


i + e> 


(1,3) 
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where (j) out = 4 > 2Lm,T2Lm is the highest momentum mode to be integrated out and the 
other two 0 are low energy modes with l < 2 L. The expectation value means that 0 out 
are contracted by use of the tree level propagator. 1/L corrections can be expressed as 
derivatives on 0, but it is sufficient to look at the leading term for the current purpose. If 
(—l) z is absent, the leading term is simply tr N (p 2 , but it cannot be organized in that way 
due to this phase factor. Inspired by the parity property of the spherical harmonics, 

Y lm {e, ip) =(-l)^ m (7T - 9, if + 7T) , (1.4) 


we may include the phase factor by the definition of a new matrix, 


—1 ) l( j>lmTlm 


(1.5) 


l.m 


and call it an antipode matrix, or an antipode held on the fuzzy sphere^ as the afore¬ 
mentioned property implies that this held resides at the antipodal point on S 2 . Thus, the 
leading term is written as tr/v(00 A ), which can be seen as the most nonlocal two-point 
interaction term on S 2 . In fact, by using (1A.8|) and (IA. 12)1 in Appendix [A] we find 


N 


trjy 


') = 


dfl 

An 


- 6. ip + it ) 


( 1 . 6 ) 


Actually, integrating out the highest modes generates various types of new nonlocal in¬ 
teraction terms with antipode helds. It is worth noticing that this nonlocality is different 
from what arises in the usual noncommutative held theory defined by the star product El 
Indeed, according to (1A. 1 5(1 in Appendix [Al a noncommutative product between func¬ 
tions 0 i (0,<p) * * 02 (6,<p) corresponds to the matrix product 0102 , which is different from 
0i (6, (p)02(7r —6 1 , tp + n) corresponding to the matrix product 0i0^. In the previous paper, 
we simply drop these new terms, and carry out RG analysis. 

However, in the hrst place, according to the spirit of RG, we have to include interac¬ 
tion terms with the antipode helds from the beginning because they are generated in the 
RG. Furthermore, as well known in quantum held theory on noncommutative geometries 
(we call it noncommutative held theory or NCFT), such an IR effect due to the UV loop 
integral is a characteristic feature of NCFT (UV/IR mixing [13]), and it is then of great 
importance and interest to investigate them in further details. If the terms with antipode 
helds are kept, the resultant effective action is not a smooth function of the momentum l 
even in the low energy regime, due to the oscillating sign. In [15], Vaidya has argued that 
this is due to UV/IR mixing effect and concluded that the Wilsonian RG cannot be imple¬ 
mented in scalar held theory on the fuzzy sphere. However, Chu, Madore, and Steinacker 
have criticized this conclusion by showing that integrating out the all momenta, instead 
of just l = 2 L, leads to much milder behavior and the two-point function is smooth at 
small values of the external momenta l [16j. They have also shown that a contribution to 
the two-point function of noncommutativity, which they call noncommutative anomaly , 

Ht is possible to define the antipode basis = (—1 ) l Ti m and write (j) A = 

*Such nonlocality may reflect a stringy nature of the theory. In a slightly different context, a string-like 
degree of freedom that connects two opposite points on a squashed fuzzy sphere has been studied in m- 


4 











remains finite even after taking a continuum S 2 limit. Thus this reflects the effect from 
the underlying noncommutativity. This noncommutative anomaly is also shown to lead 
the usual UV/IR mixing effect if we take a noncommutative plane limit by blowing up 
a point on the fuzzy sphere. In the context of the large-IV RG analysis, this observation 
suggests that we would have a well-defined RG flow by performing RG transformation 
repeatedly. If so, the large-TV RG will open up a possibility that noncommutative field 
theories are formulated constructively through large-TV matrix modelsH In this paper, we 
pursuit such possibilities and will consider the large-TV RG with antipode fields included. 


This paper is organized as follows: in the following section we formulate the large- 
TV RG based on the fuzzy spherical harmonics with antipode fields. By integrating out 
high-momentum modes perturbatively, the RG equations are derived. In Section [3l we 
find various fixed points and also carry out linearized study around them. There, we 
encounter a peculiar feature of the large-TV RG with antipode fields; namely, the RG 
equations are quite different for the cases with the starting size of TV being even or odd. 
Irrespective of this, it turns out that the position of fixed points and the critical exponents 
for each fixed points agree in even and odd TV cases. Thus, it suggests that out large-TV 
RG analysis correctly captures universal features of the theory. Finally, in Section [H we 
provide conclusions and discussions. In Appendix [V] we summarize the basics of matrix 
model formulation for scalar field theory on a fuzzy sphere and enumerate useful formulas. 
Some details of perturbative calculations are presented in Appendix [B] 

2 Large -N renormalization group on fuzzy sphere with 
antipode fields 

In this section, the formulation of the large-TV RG of the matrix model describing scalar 
field theory on the fuzzy sphere with antipode fields is presented. In the first subsection, 
we set up our action for scalar field theory on the fuzzy sphere that also includes interaction 
vertices of antipode fields. This is the starting point of our large-TV RG analysis. We then 
carry out a perturbative calculation that corresponds to a coarse-graining procedure in the 
following subsection. Finally, we explain suitable rescaling of the traces and the radius of 
the sphere pisr, with respect to various noncommutative parameters in the large-TV limit, 
and we derive a set of RG equations for mass parameters and coupling constants. 

2.1 Action with antipode fields 

In the previous work, the basic action (11.11) is used to formulate the RG equations for 
m 2 N and /% . As discussed in Section 11.11 we consider large-TV RG analysis with antipode 

II Since nonlocal interactions are introduced, apart from the noncommutative anomaly, there may be 
further points to be checked, such as unitarity or causality (in Lorentzian case, for example cm for field 
theory to be well-defined. We leave this point for future study. 
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fields. Thus, we take the following extended action as our starting point, 


C _ c( kin ') I c(P ot ') 
O J\f - D I o ]\ 


'N 
dkin.) 


J N 


Sr" =^tr N (-i[i„0] 2 + - C f\L„4,][L„^} + dfk, 

Pn 


s N ot ' ] = j^ iiN ( k !?V 4 + + «w a) 0 2 (0 A ) 2 + f\u A ) 2 '), ( 2 . 1 ) 

where m% and are mass parameters, (n a parameter for the kinetic term with an 
antipode field, and coupling constants. The index a takes the values 0,1,2a, and 
2/3 to distinguish the types of the vertices, and in particular, their numbers represent 
how many antipode fields are contained in the vertices. As shown in Appendix IA.31 
the antipode projection 0 n- c/> A has a property that all the fields inside a trace can be 
antipode projected by reversing the order of them inside the trace. See flA.241) . Thus, 
without loss of generality, we can restrict the number of the antipode fields in a trace to 
be no more than half of the number of the fields. Hence (12.11) is the most general action 
with the rotational symmetry 50(3) on the fuzzy sphere and Z 2 symmetry (j)i m — y 
up to the second derivative in quadratic terms and 0(</> 4 ) non-derivative interactions. It 
should be obvious from this observation that introducing (f) A does not mean introduction 
of new degrees of freedom, since the actual degrees of freedom are the mode coefficients 
< t>im■ This extended action is considered to be just an addition of momentum dependent 
interaction terms. On the other hand, it is much simpler to use <p A to express the action 
with such interactions since it enables us to write them in terms of the matrix product 
and trace. Thus, in the action we introduce the antipode field (j) A as if it is an independent 
field, but we note that we do not have to integrate <f> A as an independent variable of </> in 
the large -N RG. 


2.2 Coarse-graining: perturbative calculation 

In the idea of the large -N RG, the modes with high l will be integrated out to produce an 
effective action for the modes with lower angular momenta. We may start with a generic 
situation; namely, n number of the outermost momentum shells, l = 2 L, 2L — 1, • • • ,2 L — 
n + 1, are integrated out. Thus, we first divide the space of the angular momentum A 
into the one for higher modes to be integrated (called out modes), and the others (called 
in modes) as 

A = {(l, m) | 0 < l < 2L, -l < m < 1} = A f ( ”J ® A-” 0 , 

A-aui = {(l, m) | 2L — n + 1 < l < 2L, 

Aj”) = {(l, m) | 0 < l < 2L — n, —l < m < l } , (2-2) 


and define correspondingly 


*film 


<tc Me aJS 

$m (*, m ) e A in ] ’ 


(2.3) 
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A matrix only with in-modes, </> in , and similarly 0 out , are defined as 


r = E ■r* = V «T im . 

(Z,m)eA i i" ) 

The coarse-graining procedure is then formulated as 


(2.4) 




= - log 


II 


out e -S N (m 2 N ,fh 2 N ,K ( “’) 


(2.5) 


6,m)eA<,"| 


where SN-h(m 2 N _ h , m 2 N _ hl K^-h) I s an action °f (N — n) x(N — n) matrix. As in the previ¬ 
ous study [lj, the large-A" RG formalism presented here respects the rotational symmetry 
50(3), and the resultant action should also be organized with respect to 50(3) irreducible 
representations. So far, we construct a mapping from N x N matrix to (N — n)x (N — n) 
matrix. Instead of integrating out l = 2 L, 2L — 1, • • • 2L — h + 1 modes in one coarse- 
graining procedure, we may repeatedly perform the smallest n — 1 mapping n times to 
realize the same mapping. As we discuss in Section 13.2.11 truncation or approximations 
make difference between them in general, but it can be negligible to lower order calcu¬ 
lations we consider. Thus, we restrict ourselves to considering n — 1 case for the time 
being, and establish large -N RG equations. 

We stress here that our RG (12.51) is not just a mathematical problem of interest, but 
has firm physical ground. In fact, the mapping rule (1A. 1 1|) (1A. 13[) given in Appendix [A] 
tells us that the action (12.ip we start from is completely equivalent to a field theory on a 
fuzzy sphere with radius p 2 N 
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p 2 N dfl 

An 


- (£i<l >( 0 , <A>) + ^0(0, V) 2 

~ 7ri"£i</>(0, <p) A0(tt - 9, ip + n) + ^0(6 I , p)cj)(n - 9,<p + n) 
zp N z 

(0) .(1) _(2a) 

+ -^0(0, pf + pf^ - 0, p + tt) + 0(0, p) 2 ^ ~9,p + n) 2 

(2/3) 2 \ 

+ (0(0, vO0(tt - 0 , p + n )) J , (2.6) 


as a function of Here all products are understood as a noncommutative product 
defined in (1A.15I) 1**1 Hence (12. 5 j) can also be regarded as applying the standard Wilsonian 
RG to the field theory (j2,6[) on the fuzzy sphere with antipodal interactions, which is cut 
off in such a way that the rotational symmetry is preserved. Thus our RG also reveals 
properties of such a field theory. Furthermore, as we discuss later, our RG is formulated 
to describe the large-.A limit with noncommutativity a of the fuzzy sphere fixed (see (IA.5D 
in Appendix [A]) . Therefore we can fix a as small as we like to describe a field theory on 
the sphere regularized in a rotationally invariant way by small enough noncommutativity. 

**As is well known, 0((f> 2 ) or 0(4>(j ) A ) terms are not affected by the noncommutativity and hence the 
product can be replaced with the usual one as in m - 
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To carry out the integration for out modes by perturbation theory, we first decompose 
the quadratic part of the action as 


Q(kin.) _Q(kin.)in . Q(kin.)out 


s 


„2 ~ 2 


(kin.) in _ 1 ^_ ( l fr jm]2 , PN m N j.in 2 Cn r r j.ini r r j.inT , PN m N j. infill- 4 


N =xtoN\ ~o[Li, <f> m ] + 


2 

„2 


4 km ' )OUt =^tr,v ( —-[Lj, 0 out ] 2 + (p out ][Li, r utA ] + 


2 [L i: r][Li,r}+ 2 

Cn I 


„2 ~ 2 


2L 


?[ N ( N - 1 )( 1 + l -1 )^ 1 C N )+P 2 N ( m 2 N + (-l) N 1 rh 2 N )]4>^<fi. 


out 
2 L m 


m=—2L 


(2.7) 


Note that there is no cross term between 0 in and 0 out in the quadratic part of the action. 
Here, we have used the fact that for the out modes the antipode projection is simply 
multiplying the phase factor (—l)^ -1 , 


ut^ 


2 L 


E (-l) 2I 02ir.Ik™ = (-1)"-^°"* 


( 2 , 8 ) 


m=—2L 


This fact is also useful to organize the interaction vertices below. 
We define 


2 L 


Z 0 = 


FI ^ 


rvfkin.) out 

out „—S\, 


2 L 


2 Lm 


• <°)o = / n 


(-((kin.) out 


(2.9) 


m=—2L 

then our RG equation (j2.5[) becomes 
,2 ^.2 


m=—2L 


= -logZ 0 + - log(e ^ 


—5 


(pot.) 


( 2 . 10 ) 


Thus the calculation of 5 'jv-i(^at_ 1 , k n-i) amounts to evaluating (e s n ^ . Now, in 
order to carry out perturbative calculation, we reorganize the interaction part of the action 
4 pot ' } according to the number of 0 m , </> out as well as the number of antipode projections, 


o(pot-) _Pn 
bN ~N 


E «vr + E ^R’ + El^d^+^rT 

L a=0,1,2a, 2/3 6=0,la,1/3 e=0,l 

+ 3 >f , + 9 >r 


( 2 , 11 ) 


where 










M = jtr N , v< 2 ® = hr> (AM) 2 ), 


( 2 . 12 ) 

Vj 0) = rr.v (M 0 ” 1 ) , Vj'“> = itr N (>“• (<T V*" 4 + M" 2 )) , 

v{ 1/3) = trjv (0 out 0 in 0 inA 0 in ) , (2.13) 

M = tr„ (rV“ 2 ) , v< 0)wp = ltr,v (MW) , (2.14) 

M = hi\.y (MAM + .TV)) . Vf )wp = itrjv (0 O ” , O OU V‘°' 4 ) , (2.15) 

Vf = l.r'.v (o“ , o° ,ll3 j , (2.16) 

Vf = jtr w (M) , (2.17) 


and the combinations of the original couplings are packed into g^j as 


9 { on = ^n\ (a = 0,1, 2a, 2/3) 


(2.18) 


fO)_(°) . 1/ i-iAf-l (1) „(!<*)_ 1(1) , I -i \N—1 (2a) (1/3) _ 1 (1) , / , NiV-1 (2/3) 


+ 7 (-i) iV “ 1 «)v, <77 = ^7 + (-7 


ft 


TV 


9\7 = ~ a k n + (-i 


N ) 


„(°) P - JO) , W , 1 (2a) 

i/2TV — K N + 2^ " ' 2^ n ’ 

,(1)P _ 1 M) 1 ( 1\N-lj2a) 


(2.19) 

9 ( 2n NP = + ^(-i)^ -1 ^ + k% p) , (2.20) 


AT—1 ,.(2,3) 


9^=2 K n+ jt- 1 ) +(-!)" 


(l)Arp _ - li) , 

TV 1 92N — 2^ n + V " ^ 


.( 1 ) 


TV-1,d 2 ") 
'JV > 


( 2 . 21 ) 

( 2 . 22 ) 


The numbers in the lower indices for V7 and stand for the number of out modes; 
namely the number of “legs” in the following perturbation theory. On the other hand, the 
numbers in the upper indices denote the number of in fields with the antipode projection. 
The other labels, a, (3, and P and NP , are for distinction of the types of vertices. The 
symmetry factors in (l2.12jMI2.17j) are assigned by paying attention to the fact that 0 and 
( f> A are not independent as noted at the end of Section 12.11 Notice that since the antipode 
projection on out fields simply provides an overall alternating phase (—l)-^ -1 , only 0 m 
appears in the above vertices. We have also utilized the trace property of antipode fields, 
fjA.27j) . to reduce the number of antipode fields. 

The expectation value can be evaluated by using the propagator of the out modes 


where P/v does not depend on m, m' and has the form 

1 


Pn = 


N(N - 1)[1 + (Mj^Civ] +p 2 N [m 2 N + (-l)^- 1 ™^]' 


(2.23) 


(2.24) 
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By using this, we can perturbatively integrate out 0 out . It can be schematically summa¬ 
rized as 


- log ( e 


-S 


(pot.) 


o N 


a=0,1,2a,2/3 

1 { Pn 

2 \ N 


E <W + ^EE<$( V . 


(a) /, ; (a) 


2=1 a 

,(«)„(&) /-»)(“)■> ;( b ) 


E E S W (E ) E ) ) C + °^n) , (2.25) 

2,7 = 1 a,6 


where (• • -) c means taking the connected part. The summation over the indices a and b 
is understood as running over the possible values including P and NP for V ^ given in 
(12.14)1 and (12.1511 . 0(k 3 n ) stands for the third or higher order corrections in k)/ . Vq 
does not contain 0 out and is then directly inherited to N — 1 theory. 

The action (12.ip has a Z 2 symmetry 0 —* —0 (or more precisely, 0/ m — » — 0; m ), and 
so does 4 km>ut . Thus, expectation values that contain odd number of out modes vanish 
identically. Furthermore, and does not include 0 m and they do not 

affect the renormalization of the parameters. With these considerations, we have 


Sn-i 



a,b a 


+ 0(k 3 n ) + (0 in independent terms) — In Z 0 . (2.26) 

Here, — In Z$ is also a 0 m independent term and we will no longer write the last two terms 
explicitly. S™ = S'^ m ' m + ^ the original action in (12. ip with 0 replaced 

with 0 in . 

We first argue that in the low energy regime of 0 m fields, namely in which angular 
momenta l associated with all of 0 m are much smaller compared to 2 L, l <P L, correc¬ 
tions involving v\ n ^ and vE i.e. the second line in (12.26)1 . are negligible in the large- N 
(therefore large-L) limit. In order to evaluate these terms, it is sufficient to consider 

<tr JV (O 1 0 out )tr iV (O 2 0 out )> c , (tr iY (O 1 0 out )tr N (0 in 0 out3 )) c , (tr JV (0 in 0 out3 )tr N (0 in 0 out3 )) c , 

(2.27) 

where Oi (i — 1, 2) is a cubic order homogeneous polynomial in 0 in and 0 mA . For example, 
V[ ia) corresponds to the choice 0\ = |(0 in2 0 mA + 0 in ' 4 0 m2 ) as in (I2.13p . Since all three 
fields in Oi are in-modes, their angular momenta lj (j = 1,2,3) are all small compared 
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to the cutoff, lj <1. The total angular momentum of O t is, by the usual addition rule, 
bounded by Zi + Z 2 + h which is again much smaller than L. The trace trjv(0j</> out ) is 
nonvanishing only when the momenta of 0* and </> out are equal, and this condition cannot 
be met. Thus, this vertex does not contribute to the perturbative calculation we consider 
now, and the first two terms in (12.271) are indeed negligible. The third term has already 
been considered in the previous study p] as (V 2 ) c . We simply quote the result as 

^tr iV (0 i > out3 )tr A r(0 in 0 out3 )^ = (polynomial in L) x e 3Lln * , (2.28) 

which is exponentially small for large-L. Some more details are given in Appendix iBl 
Therefore, we need to consider 


S 


N -1 


_ Qin 1 Pn \ ' 

° N + ~ / . 


N 


<x=0,l L 


(a)P 

92N 


V. 


(a)P 


, (a)JVP 

+ m 


v. 


( a)NP\ 

/ C 


\(P\ 
2 1 N 


(0 )P (0 )P 

9-in 92N 


V (0)F V (0)P 


+2 sswr 


v *,o)P v (mr 


+ sTAST (V^ NF VT P ) C + 2 9« P (V<"» P V<" P ) C 

+2 S ‘r«s A ' p (vrvr p ) c +2 9 <r p 9 sr (vr p vr p ) c 

+ 2 S« wp (d 1)p vf pp )° + 9 <!> p 9 <!> p (v) 1)p v« p ) c 

, /,,(1)P V (1)IW\ , (1)»P (1)»P /,,(1)»P,,(1)»P\ 

“T z 92N 92N V K 2 k 2 / -I- i/2AT i/27V \ V 2 v 2 J 

+ 2J« (vf P V<“>) c + 2 9 <r P sS (vf Kp vP) c 

+ 2 9 «!> p aS (v« 1)p vf ) c + 2 9 « kp 9 <5 (vp^vf)^ 

+ 0(/t 3 ) + (irrelevant or negligible) . (2.29) 


In order to evaluate relevant expectation values, it is sufficient to consider the following 
pieces, 


tr* (0i 0 out 02 0° ut )) o = (2 N - 1) P N tr N 
trjv (0 1 0 out 020 out ) tr/v (0° ut4 ) 


010 2 A - ^0!(-A)0 2 a + 0(N - 2 ) 


= N(2N - iyP^tr N 


0i0 2 4 - ^0!(-A)0^ + 0(N~ 2 ) 


(2.30) 


(2.31) 


(tr* (0^02^) trjv (0 3 0 out 040 out )> 
= N(2N - l)P^tr N 1 


0 A 0 2 0 3 A 0 4 - — ( « J] A«(0 a 0 2 0 3 a 0 4 ) + 0 a 0 2 A(0 3 a 0 4 ) 
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+ (C>3 o O 4) + 0(N ~ 2 ) 


(2.32) 


where Oi (i = 1,2,3,4) are polynomials of 0 m and 0 111 " 4 or an identity 1, and some 
exponentially small terms are neglected. A is defined as in (lA.IOp and A 1 -'') acts only on 
Oi or Of. These formulas are derived in Appendix |B] (h = 1 case). By choosing suitable 
Oi , we can represent the various types of vertices. For example, 

V ( f )P =tv N (0i(j) 0Ut 0 2 (j) 0Ut ) with 0 1 = ^O in 0 inA + </> inA 0 in ), 0 2 — 1. (2.33) 

Then we can apply the general formulas to evaluate the expectation values. With these 
formulas, we find 


S N -1 


2 N 
Ptv r 
2 N 

pf 

2 N 

Pn 
4 N 

Pn_ 

4 N 


■trjy ([Li, 4> m ] 2 + C N [Li, (fi m ][L i: 4> mA ) 


ml + Bi(iV) US + bSTd (1 - plB 2 (N)gfl 


+ B,(N) (j<« p + UT P ) (1 - f>%B 2 (N)g 


trjv(0 m4 ) 


9m 


tr Ar(0 m2 ) 

tr JV (0 i > inyl ) 


4% - p 2 nb 2 (n) te p + h<r p 




v) ( + Ar p ) ( 9 sr+ksr) 


2 ' 

,(dp 1 


ti N {r*r A ) 


Pn (2a) _ 2 p> ( l\[\( ( n (0)P (1)NP\~ ( (1)P fjO)Np\ _ *- , (0)NPy 

5W Pn^2\JN ) II + 2^2iV ) + I 9 2 N + 2^27V ) ^\92N ) 

x tr 7V (0 in2 (0 inA ) 2 ) 

)&« - pib 2 (n) (Ur + Ut p ) 2 +1 < 4 r p ) 2 ' 


Pn r 


4A 




0(k 3 ) + (irrelevant or negligible), 


(2.34) 


where subleading contributions of 1/N are dropped. We have defined 

B\(N) m 2 N , fn 2 N ) — 2(2N — 1)P N , B 2 (N) = B 2 (N- m 2 N , m 2 N ) = 2{2N - l)P 2 N . 

(2.35) 


As shown here, we sometimes omit the mass dependence from B\(N) and B 2 (N ) to make 
expressions concise. Note that the mass dependence comes through the propagator factor 
P N given in (12.241) . 

The coefficients of each operator will be identified with new mass and coupling param¬ 
eters of the size N — 1 theory. However, the trace is still defined in the space ofjVxiV 
matrices, and the length scale p^ may also be renormalized in the spirit of Wilsonian RG. 
In the next subsection, we deal with them. 
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2.3 Mapping the trace and rescaling 

The result of the perturbative calculation (I2.34p is yet to be identified with a theory of 
(N — 1) x (N — 1) matrices. The trace is still defined in N x N space, and the matrix 
basis should be replaced with the one with a smaller size. Furthermore, after integrating 
out the higher momentum modes, the range of the momenta is changed, and it needs 
to be scaled to the original range as in the standard RG of field theory. This procedure 
involves the renormalization of the radius pn■ Together with these procedures, in order 
to fix the overall scale, we will normalize 0 in so that the kinetic term has the canonical 
normalization. Then we can fix relations between the parameters in the original theory 
and the renormalized theory. 

We first consider the mapping of the trace in N x N matrix space into (N — 1) x (N — l) 
one. This is essentially the same procedure we took in [TJ. Let us write the basis for N x N 
as T/jp (0 < l < 2 L) and (N — 1) x (N — 1) as ^ with l — 0, • ■ • , 2L — 1. Since 0)^ 
does not have components l = 2L, we can define 0)(^ = c0; m with 0 < l < 2L — 1, and 
also an (N — 1) x (N — 1) matrix as 0 = 4>i m T^ l> . Here, c is a constant to be fixed. 
By flA.Sp and (1A.10I) in Appendix [A] the kinetic term is mapped as 

^tr N ^ - [L^ 0 m ] 2 ^ = 1 tr JV _i ^ - [Li, 0] 2 ^ , (2.36) 


where Lj is the SU( 2) generator of spin L — | representation. We require that the kinetic 
term stays canonical through the RG procedure, and then c is fixed to be 1. The trace of 
the quadratic term also transforms simply as 

-^trjv(0 m2 ) = ^^ytr;v_i(0 2 ) • (2.37) 

It is easy to see that the same normalization applies if the trace involves antipode ones. 

On the other hand, the quartic vertices turns out to have a nontrivial factor. A trace 
with four matrices is written in terms of 3 j and 6 j symbols by (1A.19I) . A nontrivial N 
dependence comes from L = (N — l)/2 in two 6 j symbols, and we apply the following 
recursion relation from [T8]J 


a b l 
L L L 


sJ{2L + a + 1)(2L - a)(2L + b + l)(2L- b ) 


x 


2La/ (2L + l + 1)(2L — /) 


a b l 




i i + \Ja(a + 1 )b(b + 1) 


a b l 
LLL- 1 
(2.38) 


When a + b + l =even, by use of (1A.35I) and flA.36p . we find 


a b l 
LLL 


— I - 1 + — + -J— + 0(N~ 3 ' 
1 2N 8N 2 1 


L-\L-\L-\ 


(2.39) 


Here in quartic vertices a and b are the angular momenta of in modes, while l is summed 
over as in (1A.19D . When a + b + l is not even, we may roughly take the second term 
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in (12.381) to be subleading. This leads to a similar relation to (12.391) . but 0(N~ 2 ) term 
depends on a, b , and l. Thus it is natural to expect that the coefficient of 0(N~ 2 ) term 
is still 1/8, independent of a, b, and l, for this case due to the continuity, but we do not 
have a concrete expression at this moment. However, at least in the low energy regime 
a,b L, it is no doubt that the second term in (I2.38jl is of 0(1/N 2 ) compared to the 
first term. In summary, we can write down a mapping formula for quartic vertices as 


N 


= (1 + 0 (iv- 2 ; 


IV — 1 


-h- . ( r T 1 (^~ 1)/ti(JV—1)\ 

iimi 1 Izmz 1 Isms 1 (47714 / 

(2.40) 


Next, we consider the rescaling of p 2 N . In the usual Wilsonian RG of quantum field 
theory, we integrate out the momentum p for the interval A/b < p < A where A is a 
cutoff and b > 1 is a number associated with RG transformation. After integration, we 
perform a scale transformation p —>■ bp to get back to the original momentum space. It 
thus involves scale transformations of all the dimensionful quantities; especially fields get 
rescaled, and from this change we can read the scaling dimensions of the fields. We define 
a fuzzy sphere counterpart of b, called bn, as 

bl = -A-. (2.41) 

Pat-i 


So the question is how to define a scale for N — 1 theory. Since the fuzzy sphere preserves 
the rotation SO (3) symmetry, the most natural invariant is the total angular momentum 
L 2 , and the momentum squared is given by dividing it by the radius squared. Thus we 
require 


2L{2L + 1) 
Pat 


2L{2L - 1) 
Pat-i 


(2.42) 


which gives b 2 N = It should be noted that the relation between the fundamental 

scale a and the fuzzy sphere radius pN 


2 « 2 (iV 2 - 1) 

Pat ~ 1 


(2.43) 


(see Appendix is preserved by the same a, up to 0(N 2 ) corrections. Namely, this 
limit is a large -N limit with the characteristic scale of the fuzzy sphere a fixed which 
describes noncommutativity as in (lA.Sp . In this sense, we call this large -N limit the 
“fuzzy sphere limit”. Thus as stressed in (12.61) . our RG provides nonperturbative in¬ 
formation of the field theory with antipodal interactions on the fuzzy sphere with fixed 
noncommutativity a. 

As discussed in pQ, we can think of another large -N limit, which is related to 11011 - 
commutative field theory (NCFT) [HJ H] on the flat two-dimensional plane (see e.g. 
HSU2D]). Thus, we call it the NCFT limit. In this case, the large- N limit is taken with 
the noncommutativity 6 in NCFT fixed as [IB] 


A^ —> 00 with 9 


2 p| 

N 


fixed, 


(2.44) 
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which leads to b 2 N = This limit is to zoom up a tiny part of the sphere (say, the 

north pole), which will be approximated by a plane with a noncommutativity [y 1 , y 2 ] = id. 
In this case the RG is expected to describe a field theory on this plane (Moyal plane, or 
noncommutative plane). 


2.4 Renormalization group equations 


After finishing the rescaling and the mapping of the traces, (j2.34j) should be identified 
with an action of (N — 1) x (N — 1) size, up to negligible terms, 


Sn-i 


1 


IV — 1 


trjv-i 




2 , PN-l^N-l 72 CjV-1 


[L i ,4>][L i ,4> A \ + 


Pn-I^N-I 7 7A 


+ 


Pn- i 


2 2 JL 2 
K !v-d 4 + K ! v-i<i> 3 i A + K l k“\i' 2 {i A ) 2 + (4>i A ) 2 


(2.45) 


The identification of the parameters leads to the following RG equations (RGEs), 


Cn -i —C n , (2.46) 

m\_ x =b%m 2 N + b%B 1 (N)X} ll (K N ) (l - p 2 N B 2 {N)g^) + 0(ac 3 ) , (2.47) 

rh%_, =b 2 N m 2 N + }? n B 1 {N)X 2 n (k n )( 1 - p 2 N B 2 (N)g © + 0(/c 3 ), (2.48) 

=^4 - ^P^5 2 (iY)^4iv) + 0(/c 3 ), (2.49) 


with a,b = 
( (o) 

% = (/% , 


0,1, 2a, and 2/3. 

k ( 2 «) ,.( 2 / 3)5 


as 


Here we have introduced the functions of the four-vector 


v(V( K \- n (°) p i 1 n W NP -,M | 1 + 2 (~ 1 ) 

\ k n) —Q 2 n + 2^2AT — K N + ^ 

v(2)/ 1 , 1 (O)iVP 

a-AT (AC at/ —P 2 N + 2^2JV 


JV—1 


A 1 ), 


1 + (- l )^- 1 


-AC 


(2a) 

JV 


= tI « 


.( 0 ) 

V JV 


+ 


2 + (-l) 


JV-1 


_1 (O)JVP _ 1 


a^jv (Ajv) 2^2jv — 2 


AC 


2 

(o) 


-AC 


+(-ig-vr+(i+2(-i) w - i ) K r). 


JV 


+ 


-i) 


JV-l 


-AC 


(1) , t .(2/3) 


JV 


+ AC 


'JV 


43) 

JV 

^ 0) (^)=(T«( % )) 2 , 

(VjJW) =4Ai 1) (/c i v)Ai 2) (/c i v), 

at) =(t4( k Jv)) + (^AJ^Jv)) — 4 n\ k n )) , 
,(2/3)/ 5 /-42'l/ ,\2 


=(^(Kjv))"+ (d^W))"- 


(2.50) 


(2.51) 

(2.52) 

(2.53) 

(2.54) 

(2.55) 

(2.56) 


Notice that the subscript N for X$ and yff refers only to the alternating coefficients in 
their definitions. 

We note that the parameter for (j)(j) A kinetic term (n does not receive any correction. 
Thus, in the RG procedure, we can take £jv = Cn-i = ■■• = ( and ( can be arbitrary. In 
the next section, we start the fixed point analysis, and we will set ( = 0 for convenience. 
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3 Fixed point analysis 

Since the large-TV limit corresponds to performing our RG infinitely many times, it would 
be described by fixed points of the RG transformation. If they exist, for each fixed point 
we can also deduce the scaling dimensions of operators in the large -N limit from linearized 
RG transformation around the fixed point because the scaling dimension is response to 
the scale transformation. In this section, we look for fixed points of the set of the RGEs. 
If they exist, it would be strong evidence that a theory as in (12.61) with the parameters 
given by them exists consistently and nonperturbatively. Notice that for this reason the 
existence of a fixed point is striking in itself because we now allow quite nonlocal antipodal 
interactions. We then consider linearized analysis around the fixed points to determine 
the scaling dimensions. 


3.1 Fixed points for 1-step RGEs 

We first consider fixed points for the set of equations (I2.47j) - (j2.49j) . They are RGEs for the 
RG transformation from N x N theory to (.N — 1) x (.N — 1) theory, which corresponds to 
the case with n, the number of momentum shells to be integrated out, being 1. Therefore, 
we call them 1-step RGEs. 

Fixed points are obtained by setting 


m N-l 


= m N = mi 


m w-1 


= m N = m* 


„(“) _ 


k n- i — k n 


(“) _ K (a) 


(a = 0,1, 2a, 2/3 ), 
(3.1) 


and solving the relations for rn 2 , ml and ac* . As noted in the previous section, the 
parameter (n does not get renormalized and then we can consistently set (n = (n-i = 
• • ■ = (. We restrict ourselves to fixed points with ( = 0. 

Note that the mass parameters in P^, Bi(N), and B 2 (N) are all set to be fixed point 
values. We thus introduce the following notation, 


D* 
1 A 


n ~N(N- 1) + p 2 N (ml + (-1 ) N -'ml) ’ 
B{{N) =B^N- ml ml) = 2(2 N - 1 )P* N , 
B* 2 (N) =B 2 {N V mlml) = 2(2N~1)P* n 2 . 

To simplify the analysis, we use the following rescaled variables 


x 


(i) 


x 


( 2 ) 


rru 


= 




b 2 N — 1 


PnBn PnBn blp 2 N B2 ( N ) 

By use of them, the fixed point equations are written as 

X W = _ X [ l\yl 

o =vi a) ~ y£\y*) (a = 0,1, 2a, 2/3), 




b N ~ 1 ( 0 ) 
^4* 


Jr 

0 N 


(1 = 1 , 2 ), 


(3.2) 


(3,3) 
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(3.4) 

(3.5) 














where X^\y*) is given in (I2.50p with replaced with yl a> for each a and the same for 


„(°) 


(a) 


n\y*), and gfj = y 1°^ + (—1)^ 1 yi 1) + y* 2q) + yl 213 ' 1 as in (I2.22p . Since the definition of 
P^ involves ml and ml, the rescaling condition restricts the form of P ^ as 


y N 


1 — — (— 1 ) 


p 


N-lJ2) 
<+ * 


N 


N(N - 1) 


(3.6) 


Since (13.5p depends only on y[ a \ not x* \ one may find solutions to them. Then, by (13.4p . 
x^ are uniquely determined for a given set of yi a \ Together with (13.6p . the relations (13.3p 
determines the fixed points in terms of the original variables. Thus, the question boils 
down to finding solutions to (13.51) . It should also be noticed that the second term in the 
square bracket of (13.41) is 1 /N suppressed compared to the first term, 1, and then we may 
neglect that term to discuss leading order fixed points as long as gfj is of 0(1) at most. 

From the definitions of (V+ , one can see that the four equations (13. 5 p have purely 
numeric coefficients which depends on N only through an alternating sign factor (—l)^ -1 . 
Thus, solutions can be searched numerically, for N being even or odd separately, and all 
solutions will be of order 1. 

For even N, we fold the following fixed points (or lines), 


(y { ?\y£\y? a \y? p) ) =(o, 0 , 0 , 0 ), ( 0 , 0 ,-4,4), 

(t, 4 (t ± Vt), 2(t ± Vt), 2 ± 2y/t + tj , 


^ — (l + t \/21 + l) , 2 1, t, — (3 + t i y/2t + T) ^ , 

(t ± 2 Vt^l, -4 + 4 1± iy/t^T, -2 + 2 1± 2Vt^l, tj , (3.7) 


where the double sign corresponds in each solution and t is a parameter. Formally, any 
( 6 R solves the equations, but may be restricted to the range in which the values y[ a ^ 
stay real. The first one is the Gaussian fixed point, which obviously leads to ml = ml = 
kC } = 0. It is curious that the even N case allows one-parameter families of solutions, 
and then there are infinite number of fixed points irrespective of the fact that we have 
four equations of four variables. Note that the first two solutions are not included in the 
one-parameter families. 

On the other hand, for odd N, there are four fixed points found, 

(#, y ?\ yi 2a \ y {20) ) =(o, 0,0,0), (16, -64,28,20), 

/136 , 32\/2 136^2 272 , 64^2 170 , 40V2\ 

-±-, -64 +-,-±-,-±- . 

\ 9 3’ 3 ’ 9 3 ’ 9 3 J 

(3.8) 
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One can check that only the common fixed point for N even and odd cases is Gaussian 
one. So far, no parametric solution has been foundlHl 

These observations suggest that the RG flow is stable only around the Gaussian fixed 
point since we are considering the RG transformation from N x N theory to (N — 1) x 
(N — 1) one. The non-Gaussian fixed points (or lines) does not make sense since they are 
not really “fixed” along the RG transformations. The structure of fixed points are quite 
different, as N even ones include critical lines, but N odd ones not. This is in contrast 
to the case of N —> (N — 2) flow we analyze later, where N even and odd ones have the 
same number of isolated fixed points. 

Now we consider a linearized analysis around the Gaussian fixed point. Near Gaussian 
fixed point (namely, m%, fh 2 N , <C 1 and the same for the parameters of N — 1), at the 
leading order in the large -N limit we have 


2 7 2 

m N -1 — ^>n 


m 2 N + 


N 


E* 


a) K (a) 


K 


) - h 

-1 — u 


(a) 


N n N 


N 


"4-i 


= b% 


m~ N + 


N 




r\j 


N 


(3.9) 

(3.10) 


Here a = 0,1, 2cx, 2/3 and 

<S (0) =4, 5 (1) = 1+ 2(-l) 7V " 1 , S {2a) = 2 + 2(-1) n ~\ S m = 0, 

<f(°) = 2, 5 {1) = 2 + (-l) N -\ S {2a) =2(-1) n ~\ 5 m = 2 + 4(-l) 7V “ 1 . (3.11) 

The behavior of the coupling constants is trivial. We first look at critical lines of 
the mass parameters, fh(K n) and ui^Kn) which are defined by the following difference 
equations, 


2 - / ( a ) \ 7,2/2 - / { a )\\ ~2 ~ / (o) \ ,2 / ~ 2 ~ ( (n)\\ 

m N -\ ~ m[K K N _ : l) =b N (?n N - m(hrri N _ l - m{K' N _ 1 ) = b N {m N - m{K y N ’)) . 




(a) 




(3.12) 


Eqs. (13.1 Op and (I3.12p manifest that all variables have eigenvalue b 2 N of the RG trans¬ 
formation around the Gaussian fixed point. Since in the RG we have made the scale 

transformation with p —> b^p as discussed in (12.411) . this implies they all have the scaling 

dimension 2. Under the assumption that fh and m are linear functions of , these 
equations can be solved in the large- N limit for even and odd N cases separately, 

Even N\ fh — In 7V( » 4/c^ + /c^) , 

7h = In N[ - 2k$ - + 2 K^ a) + 2 , 

Odd N: fh = In N( - 4^ - 3^ ] - 4 K (2a) ) , 

fh = In 7V( — 2k$ — — 2k^ — 6k^ /3 ' > ) . (3.13) 

uAs discussed in the following subsections, the comparison between even and odd N cases should be 
done in terms of the original parameters, to* and so on. However, the structures of the set of fixed points 
are so different, and we do not expect common fixed points that make sense physically. 
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Thus, we observe that the the mass parameters have different coupling dependence near 
Gaussian fixed point for even and odd N cases. This suggests that even near Gaussian 
fixed point the critical behavior depends on evenness/oddness of N. Note that —4 In N 
coefficient of for fh is also obtained in the previous study [T], and we reproduce that 
result by setting other coupling constants to be zero, for either even or odd N. 

In order to compare the result with the one from the following 2-step study, we present 
eigenvectors of this linearized transformation in (m 2 , m 2 , n^ 2a \ basis. On 

this basis, the transformation matrix is of upper triangular form, with sextuple degenerate 
eigenvalue that corresponds to the canonical scaling dimension 2, and there are four 
eigenvectors, 

(1,0, 0,0, 0,0), (0,1, 0,0, 0,0), 

(0, 0, 5 + 4(-l)*-\ -4 - 8( 1) JV — 1 , 0, 3) , (0, 0, -1 - 2(-l) 7V " 1 , 2 - 2(-l) JV " 1 , 3, 0) , 

(3.14) 

where we have not normalized them. Unlike the cases of the following subsections, two of 
the eigenvectors are different for even and odd N. 

3.2 Fixed points for 2-step RGEs 

The observation in the previous subsection may imply that Wilsonian RG is not imple¬ 
mented on a fuzzy sphere, as Vaidya suggests in [T5j. However, as mentioned in Introduc¬ 
tion, Chu et al. |16J claim that the singular behavior due to the oscillating phase is not a 
true problem, and by carrying out loop integrals for all momenta the two point function 
does not have such an oscillating behavior. This suggests that the problem occurs since 
we consider only 1-step RG transformation, where only l = 2 L modes are integrated out. 
This motivates us to consider iterative application of RG transformation, which leads to 
theory of N x N matrices to that of (. N—n ) x ( N — h ) matrices, with n > 1. Choosing even 
n is also plausible, since evenness and oddness of N is preserved by RG transformation. 
In this subsection, we consider the first nontrivial h — 2 case, which we call 2-step RG 
transformation. Its fixed points are expected to correspond to the large- N limit taken in 
such a way that N = 2 M or N = 2 M + 1 as M —» oo. 

3.2.1 2-step RG equations 

By using the RG transformation (12.47)1 — (I2.49[) twice, we can write down 2-step RG equa¬ 
tions 

2 i2 1,2 2 

m N -2 —^N^N-l m N 

+ b 2 N b 2 n _ x ^Bi(1V; m 2 N , fh 2 N )X^\K N ) + B X {N - 1; m 2 N , mfr)x£l 1 (K N )') 

+ 0(k 2 ) , (3.15) 

~ 2 l2 l2 ~ 2 

m N -2 — 
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(3.16) 


+ b 2 N b 2 N _ i (IV; m 2 N , m 2 N )X {2] (n N ) + B 1 (N — 1; m 2 N , rh^X^ (k n )^ 

+ 0{k 2 ) , 

“N -2 ^ b N b N-l^N 

- b N b %-iPN ^-B 2 (AT; m 2 N , rh^y^( kn) + B 2 (N - 1; m 2 N , m^)^l 1 («j V )^ 

+ 0(k 3 ), (3.17) 

where we consider up to C?(k) terms for the mass parameters. When we look for a fixed 
point at the leading order in 1/N expansion, we check that all there are small and 
that 0(k 2 ) terms provide merely small corrections. In a later subsection, we carry out 
linearized analysis around fixed points, and there it is necessary to include only the leading 
corrections in the coupling constants, and then the above set of equations are sufficient 
for our purpose. Here, the mass parameters and the coupling constants on the right hand 
sides are m 2 N ,rh 2 N i and k$ . If we apply our RGE twice straightforwardly, we will have 
m2 N-i an d on the right hand sides, but they are all replaced with rn 2 N and fh 2 N 

respectively by using the RGEs again. More precisely, in the second application of the 
RGEs, there appear propagators with m| f _ 1 and m| r _ 1 , namely B\(N — l]m 2 N _ 1 ,m%_ 1 ) 
and Bo(N — 1; m| r _ 1 , m^_ 1 ). These mass parameters are to be further replaced by the 
RGEs 

m N-l = b N m N + ®( K ) , ™N-1 = bN^N + O(k) , (3.18) 

but, to the order we take now, it is sufficient to take the leading terms. This leads to 

Bi(N - =B 1 (N - 1 ■m 2 N ,m 2 N ) + 0(k) , 

Bi(N - 1; mjf, rh 2 N ) =2(2 N - 3)P N _ 1 , 

Pn " 1 = (JV _ 1)(JV - 2) + + (-1 ’ (3 ' 19) 

and B 2 (N — 1; rn 2 N , m 2 N ) = 2{2N — 3)P^ r _ 1 . Note that we have used the relation b2 N p 2 N _ x = 
p 2 N . In eqs. (13 .1 511 (13 .1 71) . xjj^( kn) (i = 1, 2) and y^^HN) mean that the N dependent 
coefficient (—l)^ 1 in (I2.50jl — (I2.56P should be replaced with (—l)^ 2 but the arguments 
are still k$. Namely, the subscript refers to the N dependence of the functions A’ J ^_ 1 (-) 

and 1 (-) themselves. Thus, on these RGEs, the mass parameter dependence appears 
through Bi coefficients, while the coupling constants dependence is from X and y. 

Before going to fixed point analysis, we take a look at the validity of 2-step RGEs just 
presented. Here, a map from N x N theory to (N — 2) x (N — 2) theory is defined by 
applying 1-step RG transformation twice, which can in principle be different from the one 
defined by integrating out both l = 2 L and l = 2L — 1 modes at the same time, if we make 
some approximations. For example, one possible difference is from a graph involving two 
out-mode propagators, where one is l = 2L mode and the other is l — 2L — 1. The latter 
does include these contributions, while the former (the ones we have just presented) does 
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not due to the low energy approximation used to derive the RG transformation. The 
reason is that in the first step RG integrating only l = 2L, all the other lines are assumed 
to have sufficiently low angular momentum l <C L by the low energy approximation, but 
in order to reproduce this graph one of them should have l — 2L — 1 in the next step 
RG, and hence the low energy approximation assumed in the first step does not work. 
Such difference can be understood as potential error terms for iterative application of RG 
transformation with low energy approximation. It is therefore important to estimate such 
difference. 

In AppendixlBl we present formulas for expectation values by integrating l = 2 L, • • • , 2 L— 
n+1 out modeso We use them to derive the RGEs from NxN theory to (N — 2) x (N — 2) 
theory by integrating out l = 2L, 2L — 1 modes. We do not present the details in this 
paper, but one can follow the calculation by use of the formulas in Appendix |B} We have 
found that, to O(k) for the masses and to 0(k 2 ) for the coupling constants, the result 
is exactly the same as fl3.15p ~ (l3.17p . The difference appears in 0(k 2 ) terms for the mass 
RGEs (I3.15P and (I3.16p . Therefore, we can trust our RGEs for h = 2 case. Interestingly, 
this observation also suggests that we may formulate h-step RGEs just by applying the 
1-step RGE n times. To the leading order in 1/N expansion, the difference can appear 
only in 0(k 2 ) terms in the mass RGEs and 0(n 3 ) terms in the coupling constant ones. 
We cannot fully justify this observation for now, as we do not have a complete formula for 
(v 2 > V 2 > '^ expectation value, but the structure of the asymptotic formula for 9 j symbols 

seems to support this conjecture. If this is the case, then we can carry out RG analysis 
to (.N — h) x (.N — n) size, where n <C N but can be very large. This completes the 
large- N RG analysis and, for example, enables us to draw the global picture of RG flows. 
However, we will leave this for future study, and now concentrate on n — 2 case. From 
the following subsection, we look for fixed points and study their properties. 


3.2.2 Fixed points 

As before, we set 


2 2 2 ~ 2 ~ 2 ~ 2 0) (a) (a) 

m N - 2 = m N = m* , m N _ 2 = m N = m*, K y N '_ 2 = «)/ = «;', 


and solve the following six equations obtained from (I3.15l) - (l3.17p . 

ml —b 2 N b 2 N _ x ml + b 2 N b 2 N _ x + B\(N - , 

ml =b 2 N b 2 N _ x ml + b 2 N b 2 N _ x + B*(N - 1 )X^_ x (k^ , 

«i a) =b 2 N b 2 N _ x K^ - b 2 N b 2 N _ lP 2 N (B* 2 {N)y^{K,) + B*(N - 1)341^* 

where 

B* X (N) —B^N] ml, ml) = 2(2 N - 1 )P* N , B*(N - 1) = 2(2 N - 3 )P* N _ X , 


(3.20) 

(3.21) 

(3.22) 

(3.23) 


HPrecisely speaking, one of the formulas is only valid for n = 2 case, due to a technical difficulty (we 
are about to mention it). But it is sufficient for the current purpose. 
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(3.24) 


B* 2 (N) =2{2N - 1 )P* N 2 , B* 2 (N - 1) = 2(2N - 3 )P* N \ , 

1 


D* 

r N 


p* 

r N-l 


N(N - 1) + p 2 N (ml + ’ 

1 


(•N -l){N-2) + p 2 N {ml - (-l)"" 1 ^) ‘ 
We introduce the following rescaled variables, 


x 


(i) 


x 


( 2 ) 


m 


mz = 


■) _ b N b N-l 1 „,(«) 


k[*' = 


* p 2 N P * N ’ * 

Py and Py^ can be written in terms of xf* (i — 1, 2) as 

1 - + (-l)^ 1 ^) 


D* 

r N-l 


N N(N - 1) 

1 1 — (xi^ + (—l) 7 ^ 1 ^) 


N(N - 1) i _ 2 (-l)^- 1 a;i 2) - |(l - (x£ 1} + (-l)*- 1 *®)) ' 
With these variables, the fixed point equations are 


(3.25) 


(3.26) 

(3.27) 


0 

0 


—x b ' 1 + Xn\u*) + 

=»!“> - yP(v.) - 


P*(1V- 1) 
P*(7V) 
B 2 *(iV-l) 
B*(N) 


X, 


(*) 


JV-l 


. v(“) 

y N—l 


(y*), 


(y*) ■ 


(3.28) 


They are again algebraic equations with numerical coefficients (for a given N). We assume 
that fixed points x b} and yi^ are of order 1, and then take the leading order coefficients 
in the large -N limit. Since -^Tn^ depends on these six equations are now coupled, 
but again we can numerically find four real solutions for even and odd N respectively. 

For even N, we hnd 



=( 0 , 0 , 0 , 0 , 0 , 0 ), 


( - 0.1671, 0.1671, 0.1902, -0.7608, 0.09582, 0.4748) , 

( - 0.06894, 0.1267, 0.08281, -0.6088, 0.1656, 0.4766) , 

( - 0.1198, -0.1270, 0.01456, 0.06177, 0.02912, 0.01820) 

(3.29) 


Let us call these four solutions e(i) (i = 0,1, 2, 3) in order. On the other hand, for odd N, 



x®yi°\yVy™y™) =( 0 , 0 , 0 , 0 , 0 , 0 ), 


( - 0.1252, 0.1252, 0.1068, -0.4274, 0.05383, 0.2667) , 

( - 0.05500, 0.1011, 0.05271, -0.3875, 0.1054, 0.3034) , 
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( - 0.1606, -0.1703, 0.02616, 0.1110, 0.05233, 0.03271) , 

(3.30) 


and we call them o(«) (i = 0,1, 2, 3). Apart from the Gaussian fixed point case, the 
numerical computation is done with 15 digits working precision, but we round the results 
to four significant figures and show them. So far, we have numerically confirmed these 
sets of solutions and it seems very unlikely that there exists another. 

At a glance, they do not seem to share common solutions, ffowever, the physical 
quantities are the original parameters m 2 N and the others. The rescaling factors in (13.251) 
also depends on N as in (12.430 . (I3.26P and thus we indeed find that the fixed point values 
of the masses and the coupling constants are the same for both e(i) and o(i) given by 

{a 2 m 2 , a 2 rhl, a 2 n//\ a 2 n^\ a 2 ^ 2a \ a 2 ^ 213 ' 1 ) 

= ( 0 , 0 , 0 , 0 , 0 , 0 ), 

( - 0.5010, 0.5010, 0.4274, -1.710, 0.2153,1.067), 

( - 0.2306, 0.4239, 0.2317, -1.703, 0.4634,1.334), 

( - 0.4826, -0.5119, 0.05909, 0.2507, 0.1182, 0.07387), (3.31) 


for i = 0,1, 2, 3 in order, ffere we have neglected subleading contributions in the 1/N 
expansion and a is the scale fixed in our RG describing the fuzzy sphere limit given in 
(I2.43p . Therefore, our 2-step RGE turns out to have the common fixed points for even and 
odd N cases, ft seems nontrivial because the fixed point equations (I3.21h (I3.23j) change 
for even and odd N. However one would anticipate agreement of the fixed points, since 
the fixed point action which describes the large-A" limit of our matrix model should not 
depend on evenness or oddness of N. Actually, with hindsight, we may have foreseen this 
agreement from the fixed point equations (13 .15 P — (13 .1 71) in the large -N limit. For example, 
we consider the leading order part of 1/N in (13. 15ft . with p 2 N , b 2 N b 2 N _ x , and the phases left 
untouched, 


N 




(l - l) 2 x b\- , )m; 


-i ) 5 


1 + $K + (- 


Xn\k*) + 




^N- l( K *) 


1 + (ml + (-1 ) N ~ 2 ml) 


On the right hand side, the phase factors in the first term are all (—l) w 1 , including 
the ones in X x \k*), while (—l)^ 2 for the second term. Thus the right hand side takes 
the same form for even and odd N. One can easily check that this is true for the other 
fixed point equations. Therefore, the fixed points in terms of the original variables should 
be the same for even and odd N. Note that this argument is more or less formal; the 
parameters have nontrivial N dependence of which we need to take care when we solve 
the fixed point equations. We have just observed that the (leading order) structure of 
the fixed point equations admits common solutions for even and odd N. At any rate, the 
fact that the even/odd N have the same fixed points strongly supports validity of our RG 
with antipode fields. 


23 























Furthermore, from the held theory point of view, our result also seems to be consistent 
with the claim in PTB]J : it is true that when we integrate only the highest mode, the result 
would not be a smooth function with different values for even and odd N as pointed out 
in [15], but integrating all modes in calculating correlation functions makes them smooth 
and completely regular in the 1/N expansion. In our large -N RG, we have only considered 
the integration over the highest modes and hence the RG itself is not smooth in the sense 
that it has the explicit oscillating phase (—l)^ -1 as found in [15] . The point here is that 
integration over all modes is realized in onr approach by looking at the fixed points. There 
the phase dependence in fact disappears and we have the well-defined fixed points. Note 
that the large -N limit, more precisely the fuzzy sphere limit we are considering is exactly 
the same as in [16] . 

We also notice that since the loop expansion parameter is in general p 2 N nffP^N as 
seen in (I2.47li - fl2.49li . the finite fixed points in (13.3111 are consistent with onr perturba¬ 
tive approach of the RG. In the following subsection, we further investigate the scaling 
dimensions of the operators associated with each fixed point, which also imply that our 
RG is legitimate enough to capture universality in the large -N limit. 


3.2.3 Linearized analysis and scaling dimensions 

In this subsection, we carry out linearized analysis around the fixed points found in the 
previous subsection. 

We start with rn 2 N RGE. By subtracting (I3.2ip from (13.15j) . and defining 5m 2 N _ 2 = 
m 2 N _ 2 — ml, Sm 2 N = m 2 N — ml and so on, we find, up to the linear order in 8m 2 N and the 
others, 


<5t71jy_2 — b N b N _ 1 8rri N 


2 1.2 


bib 


N U N -1 


^m 2 N d m 2 N + ^ Stiff d K ( a) 


x I B 1 (N- 1 miff, tn 2 N )xff (kn) + B X (N - 1; m 2 N , m 2 N )Xffff(ii N ) 


(i) 


fixed point 

(3.32) 


where m 2 N , mff , and nff on the right hand side will be replaced with their fixed point 
values after taking the derivatives. By noting that the mass derivatives only act on Bi, 
B\ and the tiff derivatives on X , we find from (13. 2 p 


Sm N _ 2 — bjybjy_^8mjy 

+ ^n-iB2(N)X^\k,){- P 2 n ) [5m 2 N + (-1 ff^Smff 
+ b 2 N b 2 N _,B*{N - 1 )4-r(^)(-Pw) [Srn 2 N - (-1 ff^Smff 


b 2 N b 2 N _ 




(a) 

N 


Bl(N)d< a) xff\^) + B*(N - l)d< a) Xffff(K*) 


(i) 


—h 2 J -) 2 Sm z 
— u N u N—l orri N 
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+ 1 [B' 2 (N)xS\y,) + B’ 2 (N - 1 )4tife.)]*«l 

+ - B 2 (N - l)4!iife.)]*?4 

+ fW-i E KWa M-hffe.) + b UN - Vd^xjiUvjWtf , (3-33) 

a 

where in the final part, we have noted that X^\k*) and (k*) are linear in and 
used the rescaled variables (13 ,25ft to express the coefficients by use of the values of the 
fixed points we have just found. It should also be noted that since the rescaling factors 
depends on m 2 and m 2 but not on n{ a \ derivative does not act on the rescaling factors. 
In the similar way, one can evaluate Siting and 8k}/_ 2 i and the result is summarized in 
the following matrix form, 

5v n _2 = {b 2 N b 2 N _ 1 1 6 + M)8v n , (3.34) 


where 8 v N = ( 8 m 2 N , 8 rh 2 N , 5 k^\ 8k§\8k^\8k%? ) ) t , and 8 v n _2 is defined in the similar 
way. 3-6 is 6 x 6 unit matrix and M is a 6 x 6 matrix whose elements are from the derivative 
part of the previous relations, 


M a 


M i2 =(-l) 


+ b*(n - i)41i (y*)}, 

[BZ(N)X%\y.) - B*(N - l)<4-i(?/*)] , 


B*(N) 
n-i i — 


B*(N) 


M ia =b 2 N b 2 N _j [B*(N)d y + B!(N - 1)9,.)4 (</*)] , 


(0 


M, 


al 


fefefet-i 1)" 

1 M-i I-s;(iv)] 2 


BUN)y%\v.) + BUN - 1)J>™ ,(!/.)) , 


(a) 


M a 2 =(-l) 


N 


-l ~{b%b 2 N _i 1)- 




B‘ 3 (N)y^ (</.) - b 3 *(^ - Wife.) 


(a) 


M, 


1 - 


AT^iV-l 


ab 


B*(N) 


BZ(N)d y myP{y.) + B* 2 {N - 


;(“) 


(3.35) 


where i = 1,2. On the left hand side, the ordered set (i,a) is understood to label 
the indices of the matrix. We have introduced £>3 (IV) = B 2 (N)P ^ and B 3 (N — 1) = 
B 2 (N — 1 )Pn-i f° r simplicity. By recalling that the N dependence of each coefficient 
factor, 


blbl-, = 1 + 0{N ~ l ), B* k (N) = B* k (N - 1) = 0{N l ~ 2k ) (k = 1, 2, 3), (3.36) 


and the fact that the values of the fixed points are of order 1, one can see that all 
the elements of M is of (9(1V _1 ). We are going to find the eigenvalues A for this RG 
transformation matrix b 2 N b 2 N _ 1 lQ + M. 

Before going, we make a remark on the eigenvalues of the above linearized transfor¬ 
mation matrix and the scaling dimensions associated with a given fixed point. Suppose 
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that there exists an operator O a of scaling dimension A at a fixed point. Near the fixed 
point, the operator receives a scale transformation through the RG transformation, 

Oa b^O A , (3.37) 

where b n = p^/P n -1 is the scale factor associated with the RG transformation defined in 
(12.41ft . As we have discussed in Section 1231 b n in general has the following form 

b 2 N =l + ^ + 0(N~ 2 ) , (3.38) 

where 7 = 2 for the fuzzy sphere limit and 7 = 1 for the NCFT limit. In the current 
consideration, we treat a sequence of the transformations, A — > N — 1 —> N — 2, and then 
the O a would have the following factor under the linearized transformation matrix, 

&X-i =1 + ^ + 0(N~ 2 ) . (3.39) 

On the other hand, since our RG transformation only changes A —> A ■**- 2 in total, the 
leading of the eigenvalue of the linearized RG should be 1. Moreover, (13.391) implies that 
the eigenvalue in general takes a form A = 1 + w/N + 0(1/N 2 ) (i.e. no fractional power 
of A) in the 1/A expansion from which we find that an operator corresponding to an 
eigenvector belonging to A has the scaling dimension A = wj 7 . More concretely, (13.341) 
gives an explicit form of the linearized RG in the 1/A expansion 

b 2 N b%_, 1 6 + M = 1 6 + + 0(N~ 2 ), (3.40) 

where 


M =2 7 1 6 + M| 0(JV _ 1} , (3.41) 

and M\ 0(n _ 1) means that we take the leading order term (namely, that of 0(1/N) as 
mentioned below (13.36ft f for each matrix element. From this form, it is manifest that 
if 0(1/N) terms are dropped, the eigenvalue of the RG transformation matrix is 1, the 
sextuple root. Thus if v is an eigenvector belonging to an eigenvalue A = l + tc/Aof the 
linearized RG in (13.401) . it follows that 


Mv — wv, (3.42) 

up to 0(1/N) corrections, and that v corresponds to an operator with scaling dimension 
A = w/'y. The general argument above predicts that w is 0( 1) and this is indeed the case 
because M is a matrix of 0( 1). Note that when the loop corrections are neglected, which 
corresponds to dropping M, the eigenvalue of M is 2y (sextuple root), and the scaling 
dimensions are all A = 2, as anticipated. 

Now, we can calculate the eigenvalue and the eigenvectors associated with each fixed 
point for A even and odd respectively. For numerical study, we also need to fix a value of 
7 from b 2 N defined above. Thus, we study 7 = 2 case (the fuzzy sphere limit) and 7 = 1 
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case (the NCFT limit) separately. It has been found that after substituting the values of 
the fixed points, the matrix elements (13.35j) of even and odd N coincides, up to very small 
numerical errors, for suitable pairs of fixed point values. Therefore, in this 2-step RGE 
case, even and odd N cases share the same scaling dimensions (critical exponents) and 
the corresponding eigenvectors (namely, associated operators) are also the same. As in 
the case of the fixed point values, we can see that this agreement comes from the structure 
of the linearized RGE. For example, from the middle expression in (13.331) . the coefficient 
of 5m 2 N (namely the matrix element Mu) reads 

-p»l -1 + bi(n - ■ 

If we consider the leading order part of 1 /N with p 2 N , b^b^^, and the phases left un¬ 
touched, it becomes 

^P 2 n^n^ ) ‘n-i _ < 4 } K ) _|__ 

N3 Lll + ^K + t-l)"" 1 *;)] 2 [1 + §i(ml + (-l)"- 2 * 2 )] 2 . 

Again, the phase factors in the first term in the square bracket are all (—l)^ -1 , including 
the ones in xjj\ k*), while (—l )^ -2 for the second term. Thus, if the fixed point values 
of k^ are the same for even and odd N, which has already been confirmed, the matrix 
element Mu is the same for even and odd N. One can check that this is also true for 
all the other matrix elements. Thus, we can conclude that even and odd N cases share 
the common eigenvalues and eigenvectors at the leading order in the large -N limit, and 
this will be true for any choice of 7 as long as fixed point values in terms of the physical 
variables coincide. 

The expression (13.35ft is suitable for numerical study. We present the eigenvalues w and 
the corresponding normalized eigenvectors on the basis {5 m 2 , 8m 2 , 8 k^°\ 5k^\ 5^ 2a \ 5 k8 2 ^) 
respectively for each fixed point e(*) (or equivalently o(*)) for * = 0,1,2,3 in order. For 
the fuzzy sphere limit (7 = 2 ), we foicfl 


c(0)(Gaussian) : 
w = 4 (sextuple) 


( 1 , 0 , 0 , 0 , 0 , 0 ), ( 0 , 1 , 0 , 0 , 0 , 0 ), 
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(0,0,1,-4, 0,3), 


VT7 


(0,0,-2, 2, 3,0), 


e(l): 

w = 4.000 (double) (0.8541, 0.1223, 0.1043, -0.4173, 0.05256, 0.2604) 

w = 2.884 ± 1.548* ( - 0.4987 ± 0.01239*, 0.4987 T 0.01239*, -0.1406, 0.5625, 

- 0.1511 =F 0.1870*, -0.2708 ± 0.1870*) 


w = 0.9568 ( - 0.6558, -0.6078, 0.3538, 0.02579, -0.1540, -0.2257) 

w = -5.308 (- 0.1434, 0.1434, 0.2011, -0.8044, 0.09024, 0.5131), 


*Here, we present numerical eigenvectors as row vectors. 
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e( 2 ): 

w = 4.000 (double) ( - 0.8572, -0.09639, -0.05270, 0.3874, -0.1054, -0.3033) 

w = 2.762 ± 1.424? (0.3712 =f 0.1930?, -0.6057, 0.08151 ± 0.002246?, -0.5079 ± 0.03613?, 

0.2629 ±0.1596?, 0.2013^0.2389?) 

w = -1.113 (- 0.1627, 0.6135, 0.2559, -0.2851, -0.4879, -0.4608) 

w = -4.687 ( - 0.06735, 0.1282, 0.1035, -0.7557, 0.1935, 0.5998), 


e(3) = 

w — 4.000 (double) (- 0.9862, 0.1435, -0.01657, -0.07029, -0.03313, -0.02071) 
w = 3.991 (0.6790, 0.7213, 0.02751, 0.1169, 0.05517, 0.03445) 

w = 3.767 (-0.8224, -0.5475, 0.01618, -0.1122, -0.09538, -0.04384) 

w = 2.667 (0.9526, -0.1993, -0.1283, 0.05703, 0.04414, 0.1767) 

w = -5.325 (0.5274, 0.5594, -0.1290, -0.5474, -0.2580, -0.1613) . (3.43) 

From the discussion above, the scaling dimensions A are given by A = w/2 for the 
fuzzy sphere limit case. Thus, the scaling dimensions associated with each fixed point are 


e(l) : A = 2.000 (double), 1.442 ± 0.7739?, 0.4784, -2.654, 

e(2) : A = 2.000 (double), 1.381 ± 0.7120?, -0.5563, -2.344, 

e(3) : A = 2.000 (double), 1.996, 1.884, 1.334, -2.662, (3.44) 


where c(0), namely Gaussian fixed point, case always has A = 2 as sextuple one, and is 
omitted from the list. From (j3.43[) and 63 .44[) . we obtain the following observations: 

• There exist doubly degenerate A = 2.000 within numerical errors for all cases. 
Numerically two eigenvectors with A = 2.000 also coincide with each other. The 
eigenvectors get significant contribution from 5m 2 N component; about 97% for e(3) 
case, and about 73% for e(l), e(2). Thus, we can naturally identify them as a 
“mass” term around each fixed point. They are the most relevant operator for all 
fixed points. 

• For e(l) and e(2), we have complex eigenvalues (thus complex scaling dimensions) 
and eigenvectors associated with them. The matrix M is a real matrix, and then 
complex eigenvalues and eigenvectors have to show up in complex conjugate pairs. 

• The e(3) fixed point has only real eigenvalues. There exist four relevant operators, 
and three of them have almost A = 2 value. It may be interesting to consider to 
define an interacting continuum field theory around it. 

We make closer look at operators with A = 2.000 and those with complex A. It is 
observed that when we estimate the fixed point values (13.3111 without enough precision, 
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two eigenvalues near 2.000 are not degenerate and eigenvectors belonging to them are 
in fact distinct. However, as the precision grows, both eigenvalues come close to 2.000 
and eigenvectors also tend to coincide, and eventually we have single eigenvector with 
doubly degenerate eigenvalue 2.000 within much precision. Physically important question 
here is whether the eigenspace with the scaling dimension 2.000 has dimension one or 
two, because it tells us the number of independence relevant operators with the scaling 
dimension of mass. The above observation seems to suggest that the eigenspace has 
dimension one and hence we have a single mass operators. At first sight, it seems strange 
that we have two types of the mass operators trjv (</> 2 ) and tr^ (< f>4 > A ) in the original 
action (12.11) . In fact, it should be noticed that given sequence of matrices for which two 
eigenvalues tend to be degenerate and their eigenvectors tend to coincide, the dimension 
of the eigenspace in the limit is not always one|H Hence we have not had a definite answer 
to this question yet for lack of the exact value of the fixed points, but we note that the 
two types of the mass term in the original action (12.ip can be written in terms of the 
modes (pi m as 


1 . (Pn™N 

tr\v 


N 




. 2 , Pn™n . 


Pn'Y] 


m 


N 


0 L + ("1 


Im'rlm 


= Pn 


E 

l : even,m 


mjy + m 


N 


0/m0? m + PN 

l : odd, m 


m N ~ m N 




(3.46) 


Thus the two mass terms can be essentially regarded as the mass terms for even l and 
odd l modes. We do not expect that they are distinct operators with scaling dimension 
of mass|j] In fact, in the fixed points in terms of physical parameters (I3.3ip as well as the 
eigenvalues and eigenvectors (13.431) . we do not find any difference between even and odd 
l and this fact ensures well-definedness of the large -N limit. Hence it may be reasonable 
to conclude that we have only one operator with A = 2 in all cases c(i) (or equivalently 
o(i)) for i = 0, • • • ,3. 

There exists a pair of complex eigenvalues in the fixed points e(l), c(2). As mentioned, 
the complex eigenvalue has to come in a pair; A = 1 + w/N with w = re™ and w* (r > 0 
and —7T < $ < 7r). By switching to a real eigenvector space, the linearized RGE matrix 
restricted to this eigenspace becomes 


pA + f 0 \ ! r /costf-sindA 

yOl + ^y 2 N ysin d cos'd ) 


tFor example, consider the sequence of matrices 



(3.47) 


( 3 . 45 ) 


with e = 1/m (m £ N). Then it is easy to see for finite m we have two distinct eigenvalues each of which 
has a one-dimensional eigenspace, and that as to —> oo, the eigenvalues and eigenvectors tend to agree. 
Nevertheless in the m —> oo limit we have a two-dimensional eigenspace of the eigenvalue 1. 

■^Since flipping the sign of rh 2 N is just exchanging the mass terms for l even and odd modes, this 
observation leads to a puzzle why there is no fixed point with the sing flipped. We do not have a definite 
answer now. 
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for a unitary matrix P. Thus, in the space of real couplings, the corresponding flow of 
RGE transformation (near the fixed point) is a spiral. When it acts on a unit vector, 
the second part gives a vector of length r/N < 1. If Kew > 0, namely cos'd > 0, this 
tiny part is outward from the fixed points, while when Kew < 0 , i.e. cos'd < 0 , it points 
toward the fixed point. Thus, the former is interpreted as a spiral source, and the latter 
is a spiral sink. When cos <p — 0 exactly, the flow of RGE transformation forms a limit 
cycle, and we need further subleading corrections in order to verify that this is the case. 
In the present case the values of d for each complex eigenvalues are 

e(l) : d = 0.4925 = 0.1568 tt , e(2) : d = 0.4760 = 0.1515 tt . (3.48) 

They are somehow close values, and correspond to spiral source behavior. 


3.2.4 Fixed points by original parameters and extra scaling factor 


So far, we have studied the properties of the fixed points for the fuzzy sphere limit case 
(7 = 2 ). However, as we have stressed in the last part of Section 13.2.21 physical RGE flow 
has to be considered in the physical parameter space spanned by m 2 N , m 2 N , k$ . When we 
look at the positions of the fixed points in this space given by (I3.25p . the finiteness of x'l ) 
and in (I3.29P and (I3.30p implies that the fixed points are located in a finite region 
in the case of the fuzzy sphere limit, where p 2 N scales like N 2 , but they go to infinity in 
the NCFT limit, N —> 00 with 6 = 2p 2 N /N fixed. The latter invalidates our perturbative 
calculation. 

In the case of well-known d = 4 — e RG analysis in scalar field theory, we need to 
include a nontrivial wave function renormalization factor for a nontrivial fixed point to be 
realized. If we canonically normalize the kinetic term, as we do here, this extra factor gives 
rise a nontrivial N dependence to the couplings. In the previous study [1], we introduced 
an extra factor c(N) = cN x attached to the coupling constant and tuned y to realize the 
NCFT limit with a finite coupling constant. Now, we try to introduce the same extra 
factor and reconsider the NCFT limit. We replace the coupling constants as 

4? c(M)4 } , (M — N — 2,N—1,N), (3.49) 


From f!3.25p . one can easily see that y should be 1 for k* to be finite in the NCFT limit. 
Then it is obvious that we cannot keep ml and ml hnite by c(lV), but as long as all 
remain hnite and small, our perturbative approach is justified. From this observation, we 
also see that we can use the common c(N ) factor for all the couplings. In RGEs fl3.15p ~ 
(I3.17p . we make the replacement (I3.49p . and do the analysis again. We now introduce a 
bit different rescaling variables, 


m 


2 

* 


ppv 


where 


~ 2 & 

m* = 


( 2 ) 


PnPn' 



£ fryfoy-i 1 ( a ) 

C (N) b 2 N b 2 N _ 1 B2(N) P H* ’ 


b 2 b 2 

U N U N -1 


c(N- 2 ) 
c(N) 


b 2 b 2 — 1 

U N U N~1 L 


7 


(3.50) 


(3.51) 
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Note that in the NCFT limit of our interest, 7 = 1 and x — 1, and then £ = 2 + 0(N^ 1 ), 
while in the fuzzy sphere limit, 7 = 2 and x — 0, and so £ = 1, where (I3.50p is reduced 
to (I3.25P (when c = 1). One can check that by use of these variables, the fixed point 
equations take the same form as (I3.28P . Thus, we can use the same set of numerical 
solutions of 0(1) given in (I3.29P and (13.30H . 

It should be noted that the propagator factor is modified too, 


P 


N 


D* 

*N-1 


1 — + (— 1 )^ 

N(N - 1) ’ 

1 1 — £(xi 1 ' ) + (—l) Ar-1 :ci 2 ' ) ) 

N(N - 1) 1 _ 2£(-l)^- 1 xi 2) - |(1 - eOri 15 + (-l)^" 1 ^)) ’ 


(3.52) 


and so are B\(N) and the other factors correspondingly. These equations with £ = 1 
reproduce (13.26p and (I3.27P . This modification does affect the position of the fixed points 
in terms of the physical parameters. By use of (13.50p . we find 


{9ml/N, 0m 2 jN , c9k^\ c0k//\ c0/4 2a) , cOk^) 

' ( - 0.4006, 0.4006, 0.1367, -0.5466, 0.06885, 0.3411) 

( - 0.1982, 0.3643, 0.08557, -0.6291, 0.1711, 0.4925) 

( - 0.4861, -0.5156, 0.02998, 0.1272, 0.05996, 0.03748) 

( - 0.5010, 0.5010, 0.2137, -0.8548, 0.1077, 0.5334) 

( - 0.2423, 0.4454, 0.1279, -0.9404, 0.2558, 0.7363) 
k ( - 0.3865, -0.4100, 0.01895, 0.08040, 0.03790,0.02369) 


e(l) solution, 
e(2) solution, 
e(3) solution, 
o(l) solution, 
o(2) solution, 
o(3) solution, 


(3.53) 


Therefore, the fixed point values do not match for even and odd N in this case. This 
may imply that this modification is not capable of defining theory around fixed points in 
a well-defined manner, but we keep moving on to the linearized analysis. At first sight it 
looks strange that the fixed point values of m 2 and rh 2 are divergent in the large- N limit 
as 0(N). However, as discussed in [T6], this is the limit that leads to a massive scalar 
held theory with the well-known phase factor associated with each vertex on the Moyal 
plane. 

The linearized RGE (I3.34p is also modified as 


5v N -2 


? I? 

’N U N -1 



+ Mj5v N , 


where M is modified from M as 


M 



(ZMij M ib \ 
\eM aj iM ab ) 



(3.54) 


(3.55) 


where i,j = 1,2 and a, b = 0,1, 2 a, 2/3, and and the others are the matrix elements 
given in (I3.35p . Note that the factors B*(N ) and B t (N — 1) are also modified as in (13.521) . 
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We again set the eigenvalue A = 1 + w/N. It is not difficult to see that, to the leading 
order in the 1/N expansion, we can drop the first and the third matrices in the definition 
of M (13.551) . Thus, we can carry out numerical study as before. It should be noted 
that we cannot forget these matrices when we fix the eigenvectors. Actually, to obtain 
well-defined eigenvectors, we need to include 1/N corrections. Thus, in the following we 
present the calculated values of the eigenvalues w, for even and odd N cases respectively. 

For even N, 

c(l) : w = 4.000, 2.044 ± 0.5627i, 2.000, 0.9568, -5.513, 

e(2) : w = 4.000, 2.618, 2.077, 1.235, -0.3694, -4.642, 

c(3) : w = 4.000, 3.470, 2.296, 2.000, 0.9338, -17.76, (3.56) 

and when N is odd, 

o(1) : w = 4.000, 2.000, 0.4404 ± 1.6837, -2.865, -16.51, 

o(2) : w = 4.000, 1.988, 0.9761 ± 1.0827, -3.180, -10.88, 

o(3) : w = 4.000, 3.769, 2.667, 2.339, 2.000, -5.673, (3.57) 

and for the common Gaussian fixed point (e(0) and o(0)), the eigenvalues are w = 2 
(double) and w = 4 (quadruple). Thus, we can observe that the scaling dimensions 
associated with fixed points do not match for even and odd N. They indeed share some 
properties; there are one A = 4 and one A ~ 2 operators for each case, for example, but 
the scaling dimensions should be universal and are expected to coincide for even and odd 
N. Together with the fact that the location of the fixed points does not agree, we conclude 
that this modification utilizing c(N) factor will not lead well-defined fixed points. 

Although it might be possible to find a well-defined NCFT limit by choosing c(N) in 
a more elaborate way, it is worth pointing out that at least perturbatively the antipode 
transformation is incompatible with the NCFT limit. Actually in this limit we have to 
restrict ourselves to the representation space of SU(2 ) with J 3 = —(N — l)/2 + 0( 1) 
corresponding to a region near the south pole of the fuzzy sphere [16]. Hence it is evident 
that even if the original matrix 0 is in this space, (ft A inevitably does not belong to 
it. In fact, we cannot figure out a possible NCFT limit of the most simple antipode 
interaction (11.61) . Since in (13.531) we found the fixed points in the NCFT limit with physical 
coupling constants as small as we can trust perturbation theory (at least by choosing c 
appropriately), this is a problem which should be addressed even perturbatively. We can 
attribute the lack of consistent fixed points and scaling dimensions to this problem. In 
[16] it is confirmed that in the NCFT limit we have been discussing, a nonplanar diagram 
on the fuzzy sphere reproduces the well-known phase factor on the Moyal plane and as 
a consequence it has IR divergence via the UV/IR mixing [14]. Since in the RG we look 
at the IR physics, it is very likely that non-existence of a well-defined fixed points in the 
NCFT limit reflects the UV/IR mixing. It is indeed true that the antipode interactions 
originate from the loop of the highest modes and we could not find any nontrivial fixed 
point including them in the NCFT limit in contrast to the fuzzy sphere limit. It would 
be interesting to examine more how the UV/IR mixing appears in our RG, in particular, 
in a nonperturbative manner. 
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4 Conclusions and discussions 


In this paper, we have formulated the large-iV renormalization group (RG) for the rank 
N matrix model which defines a scalar held theory on a fuzzy sphere. As a result of 
coarse-graining procedure in the large-iV RG, there (inevitably) appears an antipode 
held, which is defined as a scalar held twisted by a sign factor fluctuating with respect 
to its angular momentum. The antipode held is characteristic of noncommutative nature 
of the geometry. For example, it would not emerge if we regularize the theory by use of 
simple cutoff. Thus according to the spirit of the RG, we start from the action with it 
describing new nonlocal interactions between helds and antipode ones. 

It has been discussed that the appearance of the antipode held spoils the validity of 
the RG structure H3- On the other hand, it has been also shown [16] that by integrating 
momenta of intermediate states over the whole range, the renormalized action becomes 
a smooth function of the external momenta. Furthermore, such a smooth function gives 
rise to a characteristic phase factor related to UV/IR mixing effects under a suitable 
limit to noncommutative held theory. These observations lead us to expectation that 
the RG analysis provides well-defined hxed point theory that would correspond to a held 
theory on a fuzzy sphere. Actually it turns out that the renormalization group equation 
(RGE) contains the oscillating phase factor (—l)^ -1 and that the hxed points given by 
these RGEs are somehow pathological. It gives several continuum series of critical points, 
critical lines for even N , but it only provides isolated four points in the case of odd N. 
The only common hxed point is Gaussian one, and properties of Gaussian hxed point are 
also different with respect to N. Thus, we may not expect that RGE is well defined in 
this case. However, the arguments above suggest that this undesirable behavior would 
possibly be cured if we iterate the RG transformation to include the effects of integrating 
out lower momentum modes. Hence we consider an RG transformation from the original 
size N x N to (.N — 2) x (N — 2) one as a next step. This can be obtained by repetition of 
the RG transformation just given, but to the lower order in perturbation theory it agrees 
with the RGE defined by integrating out l = 2 L and 2L — 1 modes. A hxed point of the 
RGE in this case describes the large- N limit with keeping evenness or oddness of N. In 
fact, we find four hxed points for even and odd N case respectively, and also confirm that 
these four points are at the same locations. We further carry out a linearized analysis 
around them, and also observe that they provide equivalent linearized theories with the 
same scaling dimensions and the set of scaling operators. This is consistent with the claim 
in [16], because we take account of the integration over lower modes by looking at the 
hxed points and then obtain the well-defined large -N limit on the fuzzy sphere. 

More precisely, we consider two types of large -N limits; the fuzzy sphere limit in which 
the fundamental scale of the fuzzy sphere is kept, and the NCFT limit that corresponds 
to zooming up a point on a fuzzy sphere to obtain a noncommutative plane. In the case 
of the NCFT limit, we try to make the locations of the hxed points in a region where 
perturbation theory is valid by introducing an extra N dependence to couplings. This 
however does not work well; the positions of hxed points are altered differently for even 
and odd N cases, and are no longer the same. The linearized analysis also results in an 
inconsistent outcome for even and odd N. Thus, this trial may not lead to a nice NCFT 


33 


limit. This would reflect the IR singularity due the UV/IR mixing on the noncommutative 
plane, because the antipode field originates from the loop of UV modes and resulting 
interactions between fields on the antipode points are IR phenomena. Thus we find sharp 
contrast with the fuzzy sphere limit and the NCFT limit, and this observation could be 
regarded as a (nonperturbative) manifestation of the claim made in [16], that the fuzzy 
sphere does not have the UV/IR mixing, but that the noncommutative anomaly there 
yields the UV/IR mixing in the NCFT limit. 

In the usual Wilsonian RG, the more we repeat the RG transformation, the more 
kinds of interaction terms we have. In order to handle them, we usually argue that most 
of them are irrelevant and hence we could drop them. In the present case, we examine 
two-dimensional noncommutative field theories by using our large -N RG. Thus we cannot 
invoke such argument. In fact, our RG also gives rise to several derivative corrections as 
shown explicitly in e.g. (I2.30p (12.32j) . It is true that they are suppressed in 1/7V, are 
derivative terms written as the double commutator and hence can be neglected at least 
in low energy regime, but we should include them in the original action because they 
are actually generated via the RG. However, since we are in two dimensions, we have in 
principle infinitely many of them (because the scalar field has dimension zero) and it is 
impossible. Thus it is fair to say that we have found the fixed points and made analyses 
around them in the subspace of the coupling constants with such derivative corrections 
turned off. 

In our RG including antipode interactions, two well-defined fixed points have com¬ 
plex conjugate eigenvalues of the linearized RG transformation. This result is naturally 
interpreted as the fact that these fixed points are spiral sources in the two-dimensional 
subspace corresponding to operators associated with these eigenvalues. Such spiral be¬ 
havior is quite rare in the ordinary Wilsonian RG in field theories. The reasons why we 
have it are explained as follows: in order to form a spiral flow, multiple operators need 
to be mixed in the RG. However, in the ordinary Wilsonian RG we consider operators 
with definite quantum numbers like the dimensions. It is then quite hard that operators 
with different quantum numbers are mixed. In contrast, in the present case we consider 
field theories in two dimensions in which a scalar field has the vanishing scaling dimension 
at least around the Gaussian fixed point. Moreover, the most essential reason would be 
that we have exact degeneracy of operators, i.e. a field and its antipode counterpart like 
(j) and <f > A . Their degeneracy is exact and is expected to hold even nonperturbatively as 
suggested by the property given in (IA.27D in Appendix 0 Thus they would be easily 
mixed and triggers the spiral behavior. Note that scaling dimensions around a fixed point 
are usually controlled by the conformal field theory, but in our case it is not available due 
to noncommutativity. 

The other nontrivial fixed point has real scaling dimensions and eigenvectors. On 
top of that, the fixed point values of the coupling constants for the quartic terms are all 
positive]! in contrast to the other two fixed points where there exist negative ones. Thus 
it is possible that the theory around there can be well-defined, and the existence of such 
a fixed point may open a possibility that we can define an interacting field theory in a 

§Note that the values of the masses for the quadratic terms are negative, and this is quite analogous 
to the situation in usual c^ 4 theory in D < 4. 
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noncommutative space constructively. Such field theory, if any, would provide us a hint 
what kind of degrees of freedom emerges in the large -N limit from the matrices. This will 
be an interesting future direction. 

It would also be interesting to examine how nonperturbative phenomena in the large- 
N limit are captured in the large -N RG. Among others, it is quite nice if supersymmetry 
breaking shown in [22] in the matrix model describing the lower dimensional superstring 
theory [23] can be described via the large -N RG. 
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A Scalar field theory on a fuzzy sphere and useful 
formulas 

We briefly introduce scalar field theory on a fuzzy sphere, which can be represented 
as a hermitian matrix model. We also present some interesting features of antipode 
projection, introduced in Section fTTlT and also provide some formulas which are useful in 
our computation. We make the introduction concise. Readers may refer to [1] for more 
detailed introduction. 

A.l Scalar field theory on a fuzzy sphere 

We consider the following real scalar field theory on S 2 of radius p , 


s = / TT (“27 (c,mv))2 + v)2 + !' m v)i ) ' (A ' :) 


where a derivative operator £,■ = —ieijkXjdk is related to the Laplacian on a unit S 2 as 



(A.2) 


The field 0(0, <p) can be expanded by use of the spherical harmonics as 


OO 



(A.3) 


1=0 m=—l 


35 







and the action can be written in terms of the modes (j)i m . The reality condition implies that 
(j)* m = (fri-m- The spherical harmonics can be represented by use of symmetric traceless 
tensor cf m ] as 


YUO, v) = P~‘ V ■■■x\ (A.4) 

h—il 

where x l (i — 1, • • • ,3) are the standard flat coordinate of R 3 . Y,* ra = (—implies 
that = (-1 and the parity property Y lm (n - 0, ip + vr) = (-1 ) l Y lm (0,<p) is 

obvious from this expression since it corresponds to x tj ~x lj for j = 1 , • • • ,1. 

Now, we introduce a fuzzy sphere. Let L* (i = 1,2,3) be the generators of spin 
L = (N — l)/2 representation of SU( 2). We define N x N matrices x l = olL^ where a is 
a parameter of length dimension. To retain the relation x *) 2 = p 2 , a is related to p 
by p 2 = a 2 (N 2 — l)/4. Notice that 

[x\ x^] = iaeijkX k , (A.5) 

which implies that a parametrizes noncommutativity on the fuzzy sphere. Using these 
x l , we can define N x N matrices T; m , which we call the fuzzy spherical harmonics, as 

r.™ = v'E ■ (a.6) 

h’-il 

Its hermitian conjugate is T^ m = (—l) m T/ m . Let |s) be an N dimensional representation 
space of Ti m with — L < s < L. The matrix element of Tj m can be determined by 
Wigner-Eckart theorem, up to a normalization. As in [T], we use the normalization of [20] 

= (s| Ti m |s') = (-1) 1 - l j) \J (2/ + 1)N, (A.7) 


where the middle factor in the parenthesis is the Wigner’s 3j symbol. The orthogonality 
and the completeness thus follow, 


jY UiV 

S\S2^lm) S3S4 ^S1S4^- 


S2S3 1 


lm 


(A.8) 

(A.9) 


and Tim spans a complete basis for NxN matrices. Finally, corresponding to the Laplacian 
operator on S 2 (1A.2I1 . we have 


-A(j) = [Li, [Li, (j)}], where - A T Xm = 1(1 + 1 )T im ■ (A. 10 ) 

We thus define a mapping rnlc@ from a scalar held theory on S 2 to N x N matrix counter 
part as 

t Similar mapping rule for noncommutative superspace is given in m- 
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1 . function —> matrix: 


2 LI 2 LI 

EE <f>lmYlm(0, A) -t 4> = EE fafaXlm • (A. 11 ) 

1=0 m=—l 1=0 m=—l 

2 . integration —y trace: 

J ¥>) = ^ tr iv0- (A-12) 

Notice that this holds as equality. 

3. Laplacian —* double adjoint action: 

-An0(0, <p) ->■ [L i? 0]]. (A.13) 

By following these rules, scalar held theory on a fuzzy sphere that corresponds to (1A.1D 
is defined as 

s = {~F[u,<t>? + Ay + f A). (a.14) 

By putting the subscript N to p and m 2 for convenience and replacing g with k$ , we have 
our starting action m- Note that according to the mapping rule (IA. 11 j) and by using the 
fusion of T lm given in (|A.21|) later, the matrix product fa fa, with fa = J2i irni (fa)i imi Ti imi 
(i = 1 , 2 ), corresponds to a noncommutative product of fields on S 2 given as 

fa (e,<p)*fa(9,<p)= ((00 hmi L) imi l 2 m 2 (fa)l 2 m 2 ) Yl 3 m 3 (9,(p), (A.15) 

k mi (i=l,— ,3) 

where L) imi / 2m2 /l,m3 is shown in flA.211) . From this expression, we recognize that this 
star product is indeed noncommutative because the fusion Fi imi i 2m2 ll,m3 is not symmetric 
under interchange between (A,mi) and (A, ^ 2 )- 


A.2 Useful formulas of the fuzzy spherical harmonics 

In this appendix useful formulas of the N x N matrix 7] m defined in (1A.7I) are presented. 
First, the orthogonality (1A.8j) and the completeness (IA.9jl lead to 


tFv (OiT im ) tr A r { 02 ^^ — tr N (0i0 2 ), (A.16) 

lm 

— ^^trjv {oxTimP-iTi^ = IinOi tr^CA, (A.17) 

lm 


for arbitrary N x N matrices 0\, 0 2 - These properties are used to combine double traces 
into single traces. 
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By use of the formulas for the 3 j, 6 j, and 9 j symbols [18], the traces of T) m can be 
evaluated as 

tr, (1W1U1W = ** 11(2'. + 1)((-1) 21+ E(..‘. A £ M { ‘l h L H , 

tTjV iTl irni Ti^ rr i 2 Terris ) (A.18) 

4 

= A 2 J](2/ t + 1)3 ^(-l)- m (2/ + 1) 


2=1 


lm 


l\ l/± l 

mi m 


A 2 fl^h + 1)^(-1) Z2+ * 3 ^(-l)'" m (2/ + 1) 


2=1 


lm 


i h h \ lb u n hh 

m m 3 m 2 J \ L L L J \ L L L 

(A.19) 

h h l 
L L l 2 
L L l 4 
(A. 20 ) 


h h l h U 

m\ mil — m m 2 m 4 


The curly brackets with six and nine entries are the 6 j and 9 j symbols,respectively. The 
relevant formulas are also summarized in [ 1 ]. 

From (1A.8I) and flA.ISjl . we can derive the following useful fusion formula, 


rp rji _ \ ' Tp Izm^rjl 

1 hmi- L l 2 m 2 / / r l\m\l 2 m 2 1 hm 3 i 


I 3 m-3 


F 


limi l 2 m 2 


hm 3 \T~ 


JJ(2/i + l)^(-l) 2L+S: 


i =1 li+m 3 


2=1 


h h h 


m\ m 2 —m 3 / 1 L L L 


h h h 


(A. 21 ) 


By using this twice, one can easily derive a “similarity transformation” property by 


rp 

2L—n mi 
2L—n 


jLj n ( T T / 1 

X] (-1 ) m T 2L _ n m T hmi T 2L _ n _ m = A(2A -1-2 n) \ L L 2L _ n \ 


m=-2L-\-n 


{—l) n+ 1 Ti imi . 

(A. 22 ) 

If the matrix sandwiched by T 2L _ nm and T 2 ' L _ nm is an identity (namely li = 0), one finds 


2 L—n 


( —1 ) m T 2L _n m^NT 2 L-n -m —(2 A — 1 — 2?7)1jv ■ 


(A.23) 


m=—2L+n 


A.3 Properties of the antipode projection 

In this appendix, we show that the antipode projection enjoys the following property: 

(0102 •••</>m) A = ( A ‘ 24 ) 

Namely, the antipode of a string of operators is the string of the opposite ordering with 
antipode operators. 
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Proof: 

For two operators (00) A , by using the fusion formula (1A.21[) . 

(< H) A 

3 

=N* JJ(2Zi + l)^(-l) 2i+ ^=i' i+m3 
%— 1 

=^ A ct ) A . (A.25) 

Thus, 

(0102 • • • (pm) A = (02 ' • • 0™)^ fif = ' ' ' = fit ' ' ' fit fit , (A.26) 

which concludes the proof. 

From tr^r (T) m ) = NSi 0 5 m0 , it immediately follows that tr(0 A ) = tr(0). This leads to 
the following reflection property of the antipode projection inside a trace, 

trjv (01 • • • fin) =trjv (0n • • • fit) • (A.27) 

Namely, we can reverse the ordering of the fields inside a trace by putting the antipode 
projection to all of them. 


h h h 

m i m 2 —m3 


h h h 

L L L 


fih 


rni0i2m2^1 3 m 3 


A.4 Useful formulas of 3ny-symbols 

In this appendix, we summarize useful formulas we use our perturbative calculations. 
The details require more formulas than collected here, but the readers may refer to the 
previous paper p], or the textbook [18] • Many of the asymptotic relations for the 6 j 
symbols here are derived by applying the Stirling’s formula to Racah’s exact expression 
of the 6 j symbols. 


A.4.1 Asymptotic formulas and shift relations of 6 j symbols 
asymptotic formulas Racah’s asymptotic formula: for a, b , c, /, 


a b c 
baf 


cos 6 


flfl /1 (_ 1 )a+H-c + / 

b b cj ~ ^(2a + l)(2b + l) rf 
&(& ~|“ 1) ~|“ b(b -f-1) — c(c 1) 

2 \J a(a + 1 ) 6(6 + 1 ) 


P f (cos 6 ) 


Thus, for n, / -C L = (N — l)/2, 


(L L l 1 

(-ir 

[L L 2L — nj 

2 L + 1 


(-l) n [ 


N 



1(1 + l)( 2 n + 1 ) + 

2 L 

+ 1) + 0(N-‘‘) 


(A.28) 
(A.29) 


(A.30) 
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For n, m <C L, by use of the Racah formula and Stirling’s formula, 


(L L 2L — n 1 
L 2L — m J 


=2 - 4 L-2 



4 n-mrn- 2 m (~1)* m (16L 3 ) f n!m! 

[(£ + n — m)!] 2 (f!) 2 (m — £)! 


(1 + C?(i- 1 )) , 


(A.31) 


where we have assumed n > m without loss of generality due to the symmetry of the 6 j 
symbol. In the similar way, for m <C L, one can find 


2L — m 2L 2L 
L L L 


-l) 2i - m 3l(2vr)5^ f _3 /3\ 3L -^ 


8 \fm\ 


(1 + 0{L~ 1 )). (A.32) 


shift of the argument If R 1 and a, b, c are arbitrary, 
2 r{ a b c \ (-1 ) c+d+e 


(- 1 ) 


C 


c(d—e) 


d + Re + Rf + Rj ^/2R(2c+ 1) a(f ~ e) b{d ~ f) ’ 


(A.33) 


where C^ th g is the Clebsch-Gordan coefficiently whose relation to the 3 j symbol is 


CZv =(-l)“-‘ +1 v / 2?TT * _ C 7 ) . (A.34) 

When a + b + c is an even number, the Clebsch-Gordan coefficient satisfies the following 
shift property, 


a + b + c. — even : 


/^icO _ s~ycO 

^ alb-1 ~ °a-161 


c(c + 1) - a(g + 1) - b{b + 1 ) ^ e0 
2 y/a(a + 1)6(6 + 1) 


(A.35) 


By using (1X331) . one can derive the following relation 


f l h k 1 fli l 2 l 3 \ 

\lll±i j- dg l0 \lll]- 


(A.36) 


Since this is based on the asymptotic relation (1 A. 33)1 . it holds only to the leading order 
in the large-L limit. When l + l\ + I 2 is even, by (1 A. 3511 , one finds 


[ill h 1 ^ hih + 1) — l (l + 1) — h(h + 1) j l U I 2 1 
\ L L L±lj ~ 2^l{l + l)l l {l~) \ LLL ] ' 


§The formula presented in ng needs a phase factor given here. 
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A.4.2 Asymptotic formula for 9 j symbols 

We consider an asymptotic expansion formula for a 9 j symbol, 

[h h l ) 

< L L 2L — m > , (A.38) 

[L L 2 L — n J 


where l,li,l 2 ,m,n -C L. The symmetry of the 9 j symbol suggests that this is invariant 
under the simultaneous exchange of l\ yy l 2 and n m. The basic strategy to derive a 
formula is the same as p], namely by use of the decomposition into 6 j symbols, 



(A.39) 


where the triangular relations of 6 j symbols imposes the conditions, 

L — min(/, li, m) < X < L + min(7, l \), \m — n\ < l < 2L 4L — m — n , 

max(m — l, 0) < n < m + l , max(n — 1,0) < m < n + l , (A.40) 


and the usual ones for Z, Zi, and l 2 - The last two 6 j symbols are evaluated by use of the 
exact expression a la the Racah and Stirling’s formula, as presented in Appendix B in 
[lj. In the previous study, we need only n = m = 0 case. In this case, we need to take 
care of the range of the summation in Racah’s formula and the calculation is much more 
complicated. We thus work out only for the case 0 < n, m < 1, and present the result, 



(-i)"-i+V2 n h h 

\L L L 


2 N 


2 n + 1 2Zi(Zi + 1) + 2Z 2 (Z 2 + 1) — Z(Z + 1) — 1 


N 



-l) 2L+ l +1 y/2N 

16L 5 /2 


+ 1 )Iul} + 2 V /Z i( / i + 1 nLLL + 1 


+ 0(N~ 2 
(A.41) 

l 


(A.42) 


where n = 0,1. The case with m — 1 and n — 0 can be obtained by exchanging and l 2 
in the second formula. In the second formula, the shift relation (1A.37|) implies that these 
two 6 j symbols are of the same order. Thus, the second 9 j with n — 1 and m = 0 itself 
is subleading compared to the first one, namely n — m — 0,1 cases. 
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B Calculations of the expectation values 

In this appendix we evaluate generic forms of the expectation values that are necessary 
to derive renormalization group equations. We first provide general expressions for ex¬ 
pectation values by integrating out 2 L, 2L — 1, ■ ■ ■ , 2L — n + 1 out modes as in (12.211 . It is 
convenient to define a part of the out modes with the angular momentum / = 2L — n as 

= V <TV„ A-n m, (B.l) 


where n <C N. Then the out mode held (j2.3[l can be written as 0 out = J2n=o 4 > out n • Note 
that from the action (12. ip the propagator is still of diagonal form, 

(02L-nimi^2L-n 2 m 2 )o = ^n 1 n 2 ^m 1 +m 2 l) mi -P/V-ni 5 

p _ 1 

rjV-n — 


( N ~ n)(N - n - 1)[1 + (-l^-^Cw-n] + p 2 N _ n [m 2 N _ n + m 2 N _ n ] ' 

(B.2) 


We first consider 


n —1 


<tr N (Oi 4 > m 0 2 O)o = E < tr " (°i <r'J>o • 

n,m =0 


(B.3) 


where O, (i = 1, 2) are generic polynomials of (p m and . It is sufficient to calculate 
the following piece, 

(tr N (o, 02 <ry )} 0 

=S nn ,P„- n £(-l) m tr X fO, T 2 L - nm 0 2 T- 


2 L—n —m ) 


Snn'N(2N -1-2 n) P N . n E {l L 2L - n) 

l,m ^ ' 


=5 nn '(2N -1-2 n) P N - n ti N 


4 - 1 

O x O$ - P-Oi(-A)Of + 0{N ~ 2 ) 


(B.4) 


where we have used (1A.22|1 and (IA.30p . In the hnal form, O 2 can be evaluated by use of 
a property of the antipode projection (IA.24I) . — AO stands for [Li, [Lj, C?]] introduced in 

dtloli . 

Next we consider 

(tr J v(O 1 0 out O 2 0 out )tr JV (0 out ) 4 ) c 

n—1 

= E (t r »(0ir , «02r , »)tr»(r w r‘,r„r„)) I . (b.sj 

^3r" ,« 8=0 
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We need to evaluate connected graphs, which are divided into 


<tr K (O 1 ^"*„ a O 2 ^" , ni )trN(^' ,, „ s r"*„, 0 ou *„ 7 ^ ut „ s )> c 

+4i„,„ s 4„ 6 ^™(tr K (C'i«A“ ut „,C'2«A“ ut „)tr N (^"*„,<A 0 " , „4 

+4<5„,„ 7 <S„ 4 „ I <5 ws (tr K (0 1 tf> out nt 0 2 if> o '“ ni )tr N (tf> o ' lt ni tf>°'“ n5 



out 


n 5 




(B.6) 

(B.7) 

(B.8) 


where the out modes with the same index n l are to be contracted. We evaluate them 
separately. It is easy to see from the following calculation that (IB. 61) and (IB. 71) give the 
same answer, and we evaluate it as 


<tr w (O 1 ^» , n ,O 2 ^” , ni )tr N (^»* ni 0 ""‘„ s ^“*„y o ” , n5 )> o 

_ p p p \ / i \m 3 +m4+m 5 

—-'jV-n 3 -'jV-n 4 -'iV-n 5 2 _^/ V 

m 3 ,7714,777.5 

x tr n{ 0 \T 2 l —n 3 m 3 O 2 T 2 L— ri 4 7714 ) trjv (T^l — 714 — 7714 ^~2L—713 —^3 r 2L -n 5 m 3 T2L-n 3 -m 5 ) 

=<5n 3 n 4 (-PiV-n 3 ) PN-n s N 2 (2N — 1 — 2n 5 )(2N — 1 — 2n 3 ) 

A A l 


x 


E 


1,771 


L L 2L — rt 3 


-l) n+ \ 02 )l m tTN{ 0 lTi m ) 


=5 n3nA N(P N _n 3 )\2N - 1 - 2n 3 ) • P N - n5 (2N - 1 - 2n 5 ) 

A 2n 3 + 1 , a w 0 a , /nn\ T - 2 \ 


x trjy 




AT 


-C>i(-A)C>2 + 0(N~ ) 


(B.9) 


where we have again used (IA.22j) and (1A.30|h On the other hand, the rest term (IB.81) is 

(tr J v(o 1 ^" t „,o 2 ^“y)tr w (^“y^" t „y""*„ s ^" t „„)) c 

=V ,„ 4 (PN-n,) 2 PN-n s N 2 ( 2 N - 1 - 2 « 5 ) 


\m 3 


m 3 


X 


L L 2L — n 3 

A A 2A — 77,5 


(- 1 ) 2i - ri 3 + " 5t r J v(C» 1 T 2i _„3 TO3 C» 2 r 2i — n3 _ m3 ) . (B.10) 


According to the asymptotic formula (IA.31D . this is exponentially small for large-A. 
Finally, we consider 


(trjv (Oi0 out O 2 0 out ) tr^v (O 3 0 out O 4 0 out )> c 

71—1 

= E (tr K (0ir u, w 02r u y)tr K (0 3 r u, „ 7 04r u, „,)> c . (B.ll) 

«5,- ,n 8 =0 


Again, the connected pieces are 

<tr„ 2 r u *„„) tr„ (O30 ou *„ 7 O 4 0 ou ‘„,)> o 

=Ksn,K,n, (tr* (0ir u '„,0 2 r"'„ 6 ) tr„ (Oj.ryew 0 "*„J) C (B. 12 ) 
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(B.13) 


+ K„,Kn, <tr„ (Oir u *„„O 20 ”'„ e ) tr„ ( 0 3 ^"‘„ 5 04 ^ ou, „„)> c 


Due to the cyclic symmetry of the trace, the former is equal to the latter with 0 3 <-> 0 4 . 
We evaluate the latter as 

(tr„ tr N {o 3 r M H ,o,r M j) c 


=p 


N-„PN-n, J2 (-l)’” ,+mS J2 U^) 


iJknrii 


ms,me k,mi i= 1 

X trN (Xiimi^ 2 L—nsms ^27712-^21/—n6rn-6 ) (^37713^ 21 /—ns —ms^-'l^m^ 2 L—ne —me) 

4 

~-N A P N _n 5 PN-nML - 2 n 5 + 1)(4 L - 2n 6 + 1) ^ (-l) m5+me ^2/, + 1 

ms,me k,mi i= 1 

\2L— n5+/2+2L—716+/4 \ A ^ -| ^Z—m+Z'— m' ■ 


x (-1) 

'll k l 


X 


E (- 1 )' 

l,m,l' ,m' 

l 2L — 775 2L — rig 
mi m-2 mj \—fh m 5 

/1 k l ) (k k V 
x \ L L 2L — l L L 2L — n§ 

L L 2L — n§\ \lL2L — uq 


( 2 / + l)( 2 /' + l) 

h k V 


m 3 7774 m 


/' 2L — 775 2L — tiq 
-fh' —7775 —me 


~-N 4 P N _ n 5 P N _ n 6 ( 2 N - 2 n 5 - 1 )( 2 N - 2 n 6 - 1 ) e y/ 2 h + 1 

li,mi i=1 

X (_i)^+q^(_i)f-m( 2 [ + 1) / ^ k n fk k l 


1,771 

kb l ) (h k l 

x ^ L L 2L — 77.5 > < L L 2L — 77.5 

L L 2L — ns \ I L L 2L — 77-6 


mi m 2 m 1 \ m 3 m 4 — m 


(B-14) 


As noted in Appendix IA.4.21 we do not have a general asymptotic formula for this 9 j 
symbol at hand. Only available ones are those with 0 < 77 - 5 , 77-6 < 1; namely we can only 
evaluate the n — 2 case. As shown by (IA.41[) and (IA.42D . the leading order contributions 
are from 77.5 = 77-6 = 0 and 1. The result is 

(tnv {0^\0 2 ^\) tr* (0 3 r ut nO,(j>° ut n)) c 


NB 2 (N - n) 


tr^v 


o?o 2 o£o : 


2 n + 1 


2N 

+ (0 3 0 4 ) + 0(N~ 2 ) 


^AW( 0 f 0 2 0 ? 0 3 ) + 0 f 0 2 A( 0 ? 0 3 ) 

i 

, (B.15) 


where n = 0 and 1, and B 2 (N — n) = 2{2N — 2n — l) (P/v_ n ) " is used. Note that acts 
only on O t . n§ 7 ^ n§ case is 1/A suppressed compared to these contributions. 
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Finally, we cite the calculation of from (C.16) of [1], 

(tr iV ( 0 i > out 3 )tr J v( 0 i > out3 )} c 

=P 3 n (2N - 1) 3 N 4 ]T ct>? imi cj>™ m , v /(2/ 1 + l)(2/' 1 + l) 

(il,TOl),(h> m l)€Ain 


x ^(2Z + 1)(2Z' + 1) 


h 2 L l 
L L L 


l[ 2 L r 
L L L 


l 2L 2L 1 f V 2 L 2 L 1 
LL L }\L L L } 


x 


X 5 >u 


777.1+7712+777.3—m—m 7 


7771 ,••• ,7773 777 , 777 ' 


/ /1 2 L l \ 

ymi m4 m J 


( l 2L 2L\ 

V — m m3 m2 J 


x 


/ /' 2 L r \ 

y rri\ — m 2 m' y 


( V 2L 2L\ , 

\ , + (permutations of m?, 

y —??7, —m 3 —m4 y 


,?7l 4 ) 


(B.16) 


The triangular conditions for the first two 6 j symbols impose l = 2L — m and l' = 2L — n 
with m < l\ L and n < l[ L. According to the asymptotic formula of 6 j symbols 
(IA.32D . this contribution is indeed exponentially suppressed for large L, as argued in 
(I2.28p . We can therefore drop this term from our perturbative calculation. 
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